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(2)  INTRODUCTION:  Narrative  that  briefly  (one  paragraph)  describes  the  subject,  purpose 
and  scope  of  the  research. 


The  goal  of  this  research  is  to  use  high  spectral  and  spatial  resolution  (HiSS)  MR  imaging  to 
improve  images  of  human  breast.  Our  work  on  the  application  of  HiSS  to  improve  anatomic 
and  functional  imaging  was  first  described  in  a  paper  in  Academic  Radiology  [1],  Related 
work  from  this  laboratory  is  presented  in  a  number  of  other  publications  [2-9],  Work  from 
other  laboratories  shows  that  closely  related  methods  also  provide  advantages  for  anatomic 
[10]  and  functional  [11-14]  MRI. 


This  significant  body  of  work  provides  support  for  the  feasibility  of  ongoing  experiments  in 
this  laboratory.  Specifically,  we  expect  to  1)  improve  separation  of  water  and  fat  signals  2) 
increase  image  contrast  3)  increase  sensitivity  to  contrast  agents  and  to  local  physiology  - 
and  as  a  result  impiove  detection  of  suspicious  lesions  such  as  cancers  and  particularly 
delineation  of  tumor  edges.  We  expect  that  this  will  increase  the  sensitivity  and  specificity  of 

MR  scans  for  breast  cancer.  To  achieve  these  goals  our  original  'statement  of  work'  was  as 
follows: 


A.  Implementation  of  FSI  methods  on  a  clinical  Scanner:  Our  clinical  whole  body 
scanners  will  be  programmed  to  produce  oscillating  gradients  during  the  decay  of  the  proton  FID 
so  that  a  series  of  gradient  echoes  can  be  detected  following  excitation. 

B.  Processing  FSI  Data:  Spectral  information  in  FSI  data  sets  will  be  analyzed  to  reduce 
the  effects  of  resonance  offset  in  MR  images.  Then  the  corrected  FSI  data  will  be  used  to 
synthesize  images  in  which  intensity  is  proportional  to  the  peak  intensity,  linewidth,  integral,  and 
resonance  frequency  of  resonances. 

C.  A  phantom  which  contains  large  magnetic  susceptibility  gradients  and  both  lipid 
and  water  compartments  will  be  constructed  to  allow  evaluation  and  optimization  of  FSI 
methods.  Conventional  spectroscopic  images  which  use  only  phase  encoding  gradients  to  obtain 
spatial  information  will  provide  'gold  standard’  images  of  the  phantom. 

D.  Studies  of  patients:  Women  who  are  at  increased  risk  for  breast  ca  ncer  and  attend  our 
high  risk  clinic,  and  patients  who  are  treated  with  neoadjuvant  therapy  for  breast  cancer  will  be 
recruited  for  MR  studies: 

1.  Approximately  25  patients  per  year  will  be  studied  using  fast  spectroscopic 
imaging  without  contrast  agents.  FSI  will  be  correlated  quantitatively  with  conventional 
MRI  and  biopsy. 

2.  Approximately  25  patients  per  year  will  be  given  contrast.  Time  resolved  FSI 
images  of  contrast  uptake  will  be  analyzed  to  measure  rates  of  contrast  uptake  and  accurately 
identify  the  boundaries  of  enhancing  regions.  T,  -weighted  and  T2*-weighted  images  of  contrast 

agent  uptake  will  be  synthesized.  FSI  images  will  be  correlated  quantitatively  with  conventional 
MRI  and  biopsy 

3.  Quantitive  analysis  of  FSI  data  and  quantitative  comparison  with  conventional 
images:  We  will  extend  previous  work  of  Drs.  Guilhuijs  and  Giger  to  provide  quantitative 
analysis  of  FSI  and  conventional  images.  We  will  compare  edge  sharpness,  texture,  temporal 
and  spatial  gradients  in  contrast  media  uptake,  signal-to-noise  ratio,  and  contrast-to-noise  ratio  in 
FSI  and  conventional  MR  images. 

E.  MR  data  will  be  correlated  with  biopsy,  conventional  MR  images,  and 
mammography. 
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(3)  BODY: 

During  the  2  budget  year  we  made  significant  progress  towards  achieving  the  specific  aims  of 

the  proposal.  We  have  worked  towards  completion  of  all  of  the  components  of  the  statement  of 
work  listed  above: 

SOWA:  Upgraded  our  FSI  pulse  sequence  so  that  we  can  perform  fast  echo  planar 
spectroscopic  imaging  on  the  new  GE  scanners  that  were  recently  installed  at  the  University  of 
Chicago.  The  new  sequence  (referred  to  as  the  following  as  high  resolution  spectroscopic 

imaging  -  or  HiSS)  is  integrated  into  the  standard  clinical  breast  exam  so  that  it  is  more  efficient 
to  study  a  larger  number  of  patients. 

SOWB.  We  continued  development  of  methods  for  processing  the  data  and  quantitative 
comparison  of  HiSS  datasets  with  conventional  MRL 

At  the  same  time  -  we  continued  limited  studies  of  rodents  to  help  us  to  optimize  data  collection 
and  processing.  In  particular,  we  have  acquired  high  spectral  and  spatial  resolution  images  of  rat 
brain  -  because  the  well  defined  anatomy  allows  us  to  evaluate  the  data  acquisition  and 
processing.  As  before,  the  costs  of  the  rat  experiments  are  not  supported  by  this  grant,  but  the 
work  contributes  to  our  implementation  of  HiSS. 

SOWD.  We  have  increased  the  number  of  patients  and  volunteers  studies.  Specifically 
we  imaged  normal  volunteers  (n=15)  and  women  with  biopsy-confirmed  breast  cancer  (n=7) 
and  women  with  suspicious  lesions  on  mammography  (n=5).  We  demonstrated  quantitatively 
that  fat-suppression,  edge  delineation,  and  image  texture  were  improved  in  images  derived  from 
HiSS  data  compared  to  conventional  images.  HiSS  data  acquired  pre-  and  post-  contrast  media 
injection  showed  features  not  evident  in  conventional  images. 


The  work  is  described  in  detail  in  the  manuscripts  that  are  included. 


(4)  KEY  RESEARCH  ACCOMPLISHMENTS:  Bulleted  list  of  key  research 
accomplishments  emanating  from  this  research 

•  Upgrades  of  fast  HiSS  MR  imaging  methods  on  new  clinical  scanners.  The  HiSS  pulse 
sequences  can  not  be  integrated  into  standard  clinical  scans 

•  HiSS  scans  of  breast  of  healthy  volunteers,  women  with  suspicious  breast  lesions,  and 
women  who  are  being  treated  for  breast  cancer. 

•  Quantitative  analysis  of  HiSS  and  conventional  image  texture  and  edge  delineation  - 
demonstrating  that  HiSS  increases  these  measures  of  image  quality. 


(5)  REPORTABLE  OUTCOMES:: 

-  manuscripts,  abstracts,  presentations; 
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*  A  PaPer  describing  aspects  of  this  work  has  been  accepted  for  publication  in 

We  attach  the  manuscript.  This  is  an  important  step  forward  for  this  work  since  it 
will  describe  this  technology  to  Radiology’s  very  large  audience  of  practicing  academic 
Radiologists  and  hopefully  encourage  them  to  evaluate  it  in  their  own  institutions. 

fu'cAiDT1  tTted  paper  in  Academic  Radiology  (in  press)  further  demonstrates  advantages 
or  HibS  MRI  of  breast  compared  to  conventional  MRI. 


i*  T6  -are  m  t^ie  Process  °f  writing  a  second  paper  that  analyzes  the  requirement  for  high 
spectral  resolution,  and  compares  high  spectral  and  spatial  resolution  images  to  those  obtained 
using  the  Dixon  method  -  which  uses  only  two  points  of  spectral  resolution. 

•  A  third  paper  which  reports  related  research  on  the  inhomogeneous  broadening  of 
water  resonances  in  tumors  in  press  in  NMR  in  Biomedicine. 


,  T*  ^  Presented  the  work  at  last  year’s  ISMRM  and  are  submitting  two  abstracts  to  this 
year  s  ISMRM  (International  Society  of  Magnetic  Resonance  in  Medicine).  I  was  also  invited 
to  present  the  work  at  the  Contrast  Media  Research  Association  meeting  this  summer. 


List  of  papers  in  press  or  in  preparation: 

1.  Karczmar  GS,  Du  W,  Bick  U,  MacEneany  P,  Du  Y,  Fan  X,  Zamora  M,  Lipton  M; 
Spectrally  inhomogeneous  effects  of  contrast  agents  in  breast  lesions  detected  by  high 
spectral  and  spatial  resolution  MRI;  Academic  Radiology,  in  press 

2.  Du  W,  Du  Y,  Bick  U,  Fan  X,  MacEneany  P,  Zamora  M,  Medved  M,  Karczmar  GS; 
High  spectral  and  spatial  resolution  MR  imaging  of  breast  -  preliminary  experience 
Radiology,  in  press. 

3.  Medved  M,  Du  W,  Du  Y,  Bick  U,  Fan  X,  MacEneany  P.  Zamora  M,  Karczmar  G; 

Eltect  of  increased  spectral  resolution  on  water-fat  separation  in  breast  MR  imaging 
Manuscript  in  preparation  for  Journal  of  Magnetic  Resonance  Imaging  " 


-  degrees  obtained  that  are  supported  by  this  award; 

Al-Hallaq-  Ph.D  .  awarded  July,  2000 
Weiliang  Du,  PhD.  Expected  2002 


-  funding  applied  for  based  on  work  supported  by  this  award; 

1.  We  have  applied  for  a  DOD  ‘Clinical  Bridge  Award’  to  continue  and  expand 
the  work  supported  by  the  present  grant. 

2.  We  sent  a  proposal  to  the  Army’s  prostate  cancer  research  program  to  extend 
the  present  work  to  improve  imaging  of  the  prostate.  This  proposal  was  recently 
recommended  for  funding  by  the  Army  DAMD 1 7-02- 1  -0033 . 

3.  We  submitted  a  proposal  to  NIH  to  further  develop  high  spectral  and  spatial 
resolution  MRI  based  on  studies  of  rodent  tumor  models.  This  proposal  received  a  good 
priority  score  (160)  arid  it  appears  likely  that  funding  will  be  awarded. 

4.  This  work  was  a  critical  component  of  an  instrumentation  proposal  to  NIH  that 
requested  funding  for  a  whole  body  scanner.  The  scanner  would  greatly  enhance  our 
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clinical  research  capability  and  allow  us  to  scan  a  much  larger  number  of  women  with 
breast  lesions  and  greatly  reduced  cost. 


CONCLUSIONS:  Our  results  to  date  demonstrate  quantitatively  that  there  are  significant 
advantages  associated  with  high  spectral  and  spatial  resolution  imaging.  These  include 

1.  Greatly  improved  suppression  of  fat  signals  in  breast 

2.  Improved  delineation  of  edges,  for  example  tumor  boundaries 

3.  Greatly  increased  sensitivity  to  contrast  agents. 

4.  Potential  for  sensitivity  to  subvoxelar  environements  -  perhaps  microscopic  environements 
lepresented  by  the  various  components  of  inhomogeneously  broadened  water  resonances. 

During  the  coming  year  we  hope  to  expand  the  number  of  patients  and  volunteers  we  scan  and 
u  her  improve  our  methods  for  data  analysis.  We  anticipate  that  the  no.  of  patients  scanned  by 
the  end  of  the  funding  period  will  be  sufficient  to  allow  more  definitive  comparison  of  HiSS  and 
conventional  MRI  and  perhaps  a  preliminary  estimate  of  the  sensitivity  and  specificity  of  HiSS 
MRI  for  malignancies. 
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Abstract: 


Purpose:  To  evaluate  the  effect  of  varying  spectral  resolution  on  image  quality  of  high 
spectral  and  spatial  resolution  (HiSS)  imaging. 

Materials  and  Methods:  Five  women  with  suspicious  breast  lesions  and  six  healthy 
volunteers  were  scanned  using  echo-planar  spectroscopic  imaging  at  1.5  Tesla  with  0.75 
- 1  mm  in-plane  resolution,  and  2.6  Hz  spectral  resolution.  Time  domain  data  were 
truncated  to  obtain  proton  spectra  in  each  voxel  with  varying  (2.6  -  83.3  Hz)  resolution. 
Images  with  intensity  proportional  to  water  and  fat  signal  peak  heights  were  synthesized. 
Changes  in  water  signal  lineshape  following  contrast  injection  were  analyzed. 

Results:  Fat  suppression  is  optimized  at  —10  Hz  spectral  resolution,  and  is  significantly 
improved  by  removal  of  wings  of  the  fat  resonance  using  a  Lorentzian  fit.  Water 
resonance  is  often  inhomogeneously  broadened,  and  very  high  spectral  resolution  is 
necessary  to  resolve  individual  components.  High  spectral  resolution  is  required  for 
optimal  delineation  of  anatomic  features  with  very  high  T2*  (e.g.  within  a  lesion),  and  for 
detection  of  often  subtle  effects  of  contrast  agents  on  water  signal  lineshape. 

Conclusion:  Despite  a  trade-off  between  the  spectral  resolution  and  signal-to-noise  ratio, 
it  is  beneficial  to  acquire  data  at  the  highest  spectral  resolution  currently  attainable  at  1 .5 
Tesla. 

Key  Words:  high  spectral  and  spatial  resolution  spectroscopic  imaging;  echo-planar 
spectroscopic  imaging;  contrast  agent  effects;  fat  suppression;  bieast,  cancer. 


Introduction 


The  quality  of  MR  images  of  breast  has  been  enhanced  by  a  number  of  technical 
improvements  including  improved  fat  saturation  [1,2],  low  spectral  resolution 
spectroscopic  imaging  [3-11]  for  improved  fat  and  water  separation,  and  quantitative  and 
semi-quantitative  analysis  of  contrast  media  uptake  kinetics  [12-17].  Despite  these 
advances,  the  high  sensitivity  of  MRI  combined  with  limited  specificity  still  leads  to  a 
large  number  of  false  positives  [18-25],  This  is  unacceptable  if  MRI  is  to  be  used  in  a 
large  population  of  women  -  e,g.  women  who  are  at  high  risk  for  breast  cancer. 

High  spectral  and  spatial  resolution  (HiSS)  MRI  increases  the  information  content  of  MR 
data,  and  as  a  result,  has  the  potential  to  increase  sensitivity  and  specifity.  HiSS  data 
have  spatial  resolution  similar  to  that  of  conventional  MRI  or  higher  and,  in  addition,  a 
high  resolution  spectrum  of  the  water  and  fat  resonances  is  associated  with  each  image 
voxel.  Work  in  this  laboratory  demonstrated  advantages  of  HiSS  MRI  of  the  breast  [26], 
Other  work  suggests  that  HiSS  MRI  in  various  forms  -  including  multiple  gradient  echo 
imaging  [27]  -  is  advantageous  for  functional  and  anatomic  MRI  [26-33],  The 
advantages  include  improved  fat  suppression  or  fat/water  separation,  increased  anatomic 
detail  and  strong  T2*  contrast  with  lessened  artifacts  due  to  intravoxel  dephasing  effects. 
In  addition,  HiSS  provides  increased  sensitivity  to  effects  of  contrast  agents.  While 
conventional  MRI  detects  only  changes  in  T1  and/or  T2*,  HiSS  can  detect  subtle  changes 
in  water  resonance  line  shape  following  contrast  media  injection  that  may  indicate  sub- 
voxelar  regions  of  high  vascular  density.  [26, 29, 34]  Studies  of  rodent  prostate  cancer 
suggested  that  HiSS  MRI  differentiates  between  metastatic  and  non-metastatic  cancers 
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based  on  the  texture,  edges,  and  contrast  media  uptake  characteristics  in  the  primary 
tumor  [34] .  However,  high  spectral  resolution  comes  with  a  cost;  it  entails  loss  of 
signal-to-noise  ratio  and  increased  image  acquisition  time.  The  purpose  of  the  present 
work  is  to  evaluate  the  effects  of  increasing  spectral  resolution  on  various  anatomic  and 
functional  parameters  in  HiSS  images.  This  is  accomplished  by  processing  HiSS  data 
sets  at  spectral  resolutions  ranging  from  83.3  Hz  to  2.6  Hz. 

Methods 

Patient  recruitment 

Six  women  without  evidence  of  breast  abnormalities  and  5  women  with 
mammographically  suspicious  breast  lesions  were  scanned.  The  patients  were  recruited 
from  the  University  Cancer  Risk  Clinic  and  the  mammography  service.  The  reference 
(healthy)  volunteers  were  recruited  from  the  clinic  staff.  All  patients  were  studies  under 
a  protocol  approved  by  the  University  Institutional  Review  Board  after  informed  consent 
was  obtained. 

Image  acquisition 

MR  images  were  acquired  on  1.5  Tesla  GE  SIGN  A  scanners  equipped  with  ECHO 
SPEED™  GRADIENTS  (GE  Medical  Systems,  Waukeshaw  Wisconsin),  The  body  coil 
was  used  for  excitation  and  a  phased  array  coil  designed  for  breast  imaging  was  used  to 
detect  signal.  HiSS  images  were  acquired  using  echo-planar  spectroscopic  imaging 
(EPSI).  [35-37]  The  EPSI  pulse  sequence  was  composed  of  a  slice-selective  excitation 
pulse  (4  mm  slice  thickness),  a  phase  encoding  gradient  with  256  phase  encoding  steps 
and  a  train  of  128  gradient  echoes  with  256  or  384  sampled  points  per  gradient  echo. 
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This  produced  images  with  256  or  384  x  256  points  of  spatial  resolution,  with  resulting 
in-plane  resolution  of  750  -  1000  microns,  and  128  bins  of  spectral  resolution.  The 
gradient  echo  train  was  sampled  at  a  bandwidth  of  62.5  kHz,  so  that  the  FID  was  sampled 
for  approximately  400  ms.  As  a  result,  maximum  spectral  resolution  was  2.6  Hz,  or 
0.04 lppm.  The  acquisition  time  per  slice  under  these  conditions  was  approximately  130 
seconds.  The  fat  and  water  peaks  were  separated  by  more  than  64  spectral  bins,  so  that 
N/2  ghost  peaks  would  not  overlap  fat  or  water  signal.  These  ghost  peaks  are  a 
consequence  of  the  uneven  sampling  of  the  FID,  analogous  to  the  N/2  ghosting  inherent 
to  EPI  sequences.  Typically,  only  a  single  slice  was  imaged  in  each  patient.  The  patients 
were  imaged  both  pre-  and  post-  contrast  agent  injection  (0. 1  mM  per  kg  of  patient 
weight,  Omniscan™,  Nycomed  Amersham),  while  the  healthy  volunteers  did  not  receive 
contrast  agent. 

Image  processing 

The  EPSI  data  were  processed  with  a  3D  Fourier  transform  to  obtain  a  high-resolution 
proton  spectrum  in  every  voxel  of  the  imaged  slice.  An  algorithm  was  used  on  the 
resulting  data  to  detect  the  fat  and  water  peak  frequencies,  in  the  following  manner: 

1 .  The  highest  peak  in  each  voxel  was  identified,  and  a  map  of  its  frequency  was 
generated. 

2.  The  voxel  of  highest  peak  height  intensity  was  a  priori  identified  by  the  software 
as  either  fat  or  water,  and  fat  and  water  frequencies  were  determined  for  this 
voxel.  Alternatively,  one  or  more  seed  voxels  were  chosen  in  which  the  fat  and/or 
water  peaks  were  manually  identified. 
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3.  The  voxel  to  be  processed  next,  Vn,  was  selected  as  the  highest  intensity  voxel 
adjacent  to  the  already  processed  voxels. 

4.  The  frequency,  fvN,  of  the  highest  peak  in  Vn  was  compared  to  the  expected  fat 
frequency,  fexp,  equal  to  the  average  of  fat  frequencies  in  neighboring,  already 
processed  voxels.  Based  on  the  separation  of  fvN  and  feXp,  the  highest  peak  in 
VN  was  identified  as  either  fat  or  water,  and  hence  the  fat  frequency  in  Vn  was 
determined. 

5.  The  algorithm  resumed  at  point  3,  until  all  the  voxels  were  processed,  yielding  a 
map  of  the  fat  frequency  or,  equivalently,  the  Bo  map. 

In  most  cases,  this  algorithm  did  not  require  human  intervention.  It  also  effectively 
avoided  errors  due  to  Nyquist  ghosting,  inherent  in  EPI,  and  hence  in  EPSI.  In  our  EPSI 
images,  these  show  as  spectral  ghosts  due  to  uneven  sampling  of  the  FED  at  high  kx 
values. 

After  the  fat  frequency  map,  and  hence  the  water  frequency  map,  was  determined,  water 
peak  height  and  peak  integral  images  were  constructed,  as  well  as  the  analogous  fat 
images.  The  high-resolution  spectrum  in  each  voxel  was  fit  using  a  Lorentzian 
approximation  for  the  fat  and  water  peaks,  and  for  any  N/2  spectral  ghosts.  The 
Lorentzian  fits  for  fat  peaks  and  N/2  ghost  peaks  were  then  subtracted  to  allow  accurate 
analysis  of  the  water  signal.  This  was  particularly  important  for  detection  of  small 
amounts  of  water  in  the  presence  of  large  fat  signals.  As  the  water  signal  was  often 
inhomogenously  broadened,  the  peak  and  integral  water  images  were  not  obtained  from 
the  fit  parameters.  Rather,  after  subtraction  of  the  fits  for  fat  and  N/2  ghosts,  the 
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maximum  of  the  isolated  water  signal  in  each  voxel  was  used  to  construct  water  peak 
height  images,  and  the  integral  in  a  25  Hz  range  around  the  water  peak  frequency  was 
used  to  construct  integral  images. 

This  procedure  was  repeated  to  construct  water  and  fat  peak  height  images  at  varying 
spectral  resolutions  (See  Fig.  3).  The  spectral  resolution  was  varied  between  2.6  -  83.3 
Hz  by  either  using  the  full  FID  with  all  128  echoes,  or  truncating  the  FID  at  N  =  64, 32, 

16,  8  or  4  echoes,  and  substituting  the  remaining  echoes  with  zeros.  This  resulted  in 
effective  spectral  resolutions  of  2.6  Hz  (at  full  resolution),  5.2  Hz,  10.4  Hz,  20.8  Hz, 
41.6Hz  and  83.3Hz.  Truncation  was  preferred  over  exponential  or  other  smoothing,  as 
the  purpose  was  to  explore  faster  data  acquisition  with  a  reduced  number  of  echoes, 
rather  then  image  processing  methods. 

Image  characterization 

The  quality  of  the  water  peak  height  images  was  quantified  based  on  fat  suppression  and 
contrast.  The  effectiveness  of  fat/water  separation  is  assessed  by  calculating  a  fat 
suppression  ratio,  FSR.  We  define  FSR  as  the  ratio  of  the  average  intensity  of  ‘water 
peak  height  images’  in  the  voxels  containing  predominantly  water  to  the  average  signal  m 
the  voxels  containing  predominantly  fat.  [26]  Thanks  to  high  spatial  resolution  of  our 
images,  all  but  few  voxels  fall  into  one  of  these  two  categories.  Since  ideally,  only  the 
water  signal  contributes  to  these  images,  signal  in  voxels  that  contain  predominantly  fat 
should  be  minimized,  and  signal  in  voxels  that  contain  predominantly  water  should  be 
maximized.  Therefore,  effective  water/fat  separation  should  increase  FSR. 
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The  contrast  within  the  water  image  was  assessed  via  a  ‘texture’  parameter,  Sr,  which 
describes  the  non-uniformity  of  the  signal.  Sr  is  defined  as  the  increase  in  the  suiface 
area  of  measured  signal,  (when  intensity  is  plotted  against  the  2D  xy  plane  of  imaged 
slice),  relative  to  its  projection,  Sxy,  onto  this  xy  plane:  [26, 34] 

'jLifi  +  fx  +fy~S*y  (1) 

dr  -  -  * 

If  image  intensity  in  the  region  of  interest  is  constant,  Sr  is  0,  and  highly  variable  image 
intensity  results  in  a  large  5).  Only  the  voxels  within  the  breast  aie  included  in 
calculations  of  FSR  and  Sr. 

Contrast  agent  effects 

To  examine  the  effect  of  contrast  agent,  the  post-contrast  images  were  corrected  for  any 
motion  that  might  have  occurred  between  the  pre-and  post  contrast  acquisitions.  Ten 
anatomical  points  of  reference  were  selected,  and  a  least-squares  fit  was  performed  for 
parameters  of  a  linear  transformation  between  the  reference  point  coordinates  in  the  two 
images.  This  procedure  involved  some  smoothing  between  voxels,  but  dramatically 
improved  the  overlap  of  the  pre-  and  post-contrast  images.  After  the  motion  correction, 
the  pre-  and  post  contrast  water  peaks  were  corrected  for  global  B0  shifts  caused  by  the 
contrast  agent. 

The  difference  of  the  post-  and  pre-contrast  water  resonance  was  calculated  m  each 
voxel.  Voxels  with  inhomogeneously  broadened  water  resonances  were  identified  as  the 
ones  in  which  the  difference  spectra  (spectrum  post-contrast  minus  spectrum  pre- 
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contrast)  showed  a  strong  peak  whose  width  at  the  base  was  less  than  the  average  width 
of  the  pre-  and  post-contrast  peaks,  [38]  In  many  cases,  the  width  at  the  base  could  not 
be  determined,  and  was  approximated  as  twice  the  width  at  half-height.  Only  pixels  in 
which  the  pre-  and  post-contrast  resonances  were  above  the  99%  percentile  noise  level, 
and  for  which  the  difference  peak  was  above  1 ,41  x  this  value,  were  Included, 

Results 

Figure  la  shows  the  dependence  of  FSR  in  water  peak  height  images  on  number  N  of 
processed  echoes  in  the  acquired  echo  train,  for  five  patients  with  breast  lesions.  The 
frequency  resolution  increases  with  N,  as  explained  in  the  ‘Image  Processing  section. 
The  values  shown  in  diamonds  were  determined  by  fitting  the  fat  peak  with  a  Lorentzian 
function  and  subtracting  this  Lorentzian  fit  from  the  proton  spectrum  (see  Fig.  2),  to 
obtain  an  optimal  value  for  water  peak  height  in  each  voxel.  The  values  shown  in  crosses 
were  obtained  from  the  peak  height  of  the  water  resonance  without  subtracting  the 
Lorentzian  fit  of  the  fat  resonance.  For  all  patients,  the  FSR  increases  sharply  between  N 
=  4  and  N  =  32.  The  FSR  levels  off  at  N=32  (10.4  Hz  resolution)  for  all  but  one 
patient,  whose  FSR  peaks  at  N  =  16.  At  values  of  N  higher  then  32,  there  is  a  slight 
decrease  in  FSR,  probably  due  to  the  decrease  in  signal-to-noise  ratio.  Subtracting  the 
Lorentzian  approximation  to  the  fat  peak  from  the  proton  spectrum  significantly 
improves  the  FSR. 

The  advantages  of  peak  fitting  are  greatest  in  the  case  where  a  small  water  resonance 
overlaps  the  wings  of  a  very  strong  fat  resonance.  In  Fig,  2,  this  effect  is  even  more 
significant  as,  due  to  the  spectral  wraparound,  the  water  resonance  appears  at  5,2  ppm, 
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very  close  to  the  fat  peak.  (The  true  frequency  of  the  water  resonance  is  arbitrarily  set  to 
Oppm.)  The  removal  of  the  fat  signal  by  peak  fitting  allows  a  much  more  accurate 
measurement  of  the  amplitude  of  the  small  water  signal.  Off-resonance  fat  components 
are  also  evident  at  ~4ppm. 

In  Fig.  lb,  we  plot  the  dependence  of  Sr  on  N  for  the  five  patients  with  breast  lesions. 
Similar  to  Fig.  la,  we  show  values  obtained  with  (diamonds),  and  without  (crosses) 
subtraction  of  the  Lorentzian  fit  to  the  fat  resonance.  Fitting  the  fat  peak  improves  Sr 
slightly,  but  not  as  significantly  as  it  does  the  FSR.  hi  all  cases,  Sr  increases 
dramatically  with  increased  frequency  resolution,  and  levels  off  at  N  =  64  (5,2  Hz 
resolution).  However,  the  calculations  of  Sr  performed  on  the  whole  breast  may  obscure 
improvements  in  the  texture  in  smaller  but  critical  regions  with  long  T2*,  e.g,  within  the 
lesion.  When  Sr  was  calculated  over  smaller  areas  within  the  breast  lesion,  there  were 
often  improvements  beyond  N  =  64,  increasing  by  as  much  as  29%  between  N  =  64 
and  N  =  128. 

Figure  3  shows  the  water  and  fat  peak  height  images  for  a  patient  with  infiltrating  ductal 
carcinoma,  calculated  at  three  different  spectral  resolutions:  2.6  Hz  (N  =  128),  10.4  Hz 
(N  =  32)  and,  41.6  Hz  (N  =  8).  At  N  =  128,  inhomogeneities  in  the  lesion  are  clearly 
visible,  due  to  strong  T2*  contrast  in  regions  where  T2*  is  long.  This  contrast  is  reduced 
with  reduced  spectral  resolution,  which  is  particularly  apparent  at  N  =  8,  where  the 
lesion  appears  uniform. 
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Figure  4  shows  spectra  from  five  different  voxels  from  one  patient  with  infiltrating  DCIS 
and  one  healthy  volunteer  at  full,  N  =  128  (a)  and  two  lower,  N  =  64  (b)  and  32  (c) 
spectral  resolutions.  With  reduced  spectral  resolution  (b),  the  signal-to-noise  ratio  (S/N) 
is  significantly  improved,  as  would  be  expected  due  to  elimination  of  low  S/N  points  at 
the  end  of  the  FID.  However,  details  of  the  water  resonance  are  obscured.  For  example, 
at  N  =  8  (41.6  Hz  spectral  resolution)  the  water  resonance  is  significantly  broadened,  and 
truncation  artefacts  are  present.  Even  at  N  =  32  (10.4  Hz  spectral  resolution),  the  small 
shoulders  on  the  water  resonance  (arrows)  are  not  resolved.  This  may  become  an 
important  consideration  in  contrast  enhanced  MRI  (see  below).  In  pixel  #2,  the  water 
resonance  is  clearly  resolved  only  at  the  highest  spectral  resolution. 


In  Figure  5,  the  pre-  (solid  line)  and  post-contrast  spectra  (dotted  line)  from  a  single  voxel 
are  shown  at  two  different  resolutions,  at  N  =  128  (a)  and  N  =  32  (b).  At  full  spectral 
resolution,  the  inhomogenous  broadening  of  the  water  line  is  clearly  visible.  A  shoulder 
on  the  high-field  side  of  the  water  resonance  is  much  more  clearly  resolved  after  contrast 
injection  (arrow).  At  the  lower  resolution,  there  is  significant  loss  of  detail  in  the  water 
and  fat  signals,  and  the  effect  of  the  contrast  agent  is  not  accurately  represented. 

Figure  6  shows  images  that  highlight  voxels  in  which  there  are  spectrally  inhomogeneous 
changes  caused  by  contrast  agents.  The  images  are  calculated  at  3  different  spectral 
resolutions,  at  N  =  128  (a),  64  (b)  and  32  (c).  The  inhomogeneously  broadened  pixels 
tend  to  appear  near  the  boundary  of  regions  with  high  water  signals:  they  surround  the 
lesion,  and  appear  along  the  edges  of  blood  vessels.  The  number  of  highlighted  voxels  is 


largest  in  the  image  calculated  at  N  =  128  -  it  is  707  vs  593  and  451  in  images 
calculated  at  N  =  64  and  32,  respectively.  However,  the  random  placement  of  pixels 
that  emerge  at  full  spectral  resolution  indicates  that  some  of  this  increase  may  be  due  to 
the  increased  noise  level. 

Discussion 

Our  results  suggest  that  in  human  breast  at  1 .5  Tesla,  2,6  Hz  spectral  resolution  or  better 
is  required  to  resolve  distinct  components  of  inhomogeneously  broadened  water 
resonances  (Fig,  4).  In  addition,  very  narrow  water  lines  can  be  accurately  measured, 
providing  T2*  contrast  within  a  lesion,  which  is  obscured  at  lower  spectral  resolutions 
(Fig.  3).  The  intra-lesion  ihnomogeneities  are  important  for  lesion  characterization,  and 
thus  the  obtained  T2*  contrast  is  diagnostically  valuable. 

Spectrally  inhomogeneous  effects  of  contrast  agents  are  resolvable  at  resolutions  of  10.4 
Hz  or  better,  allowing  a  novel  approach  to  the  use  of  contrast  agents.  While  conventional 

MRI  measures  only  T1  and  T2*  changes  following  contrast  agent  injection,  it  may  be 
useful  to  measure  changes  in  water  signal  lineshape.  These  effects  may  be  only  weakly 
connected  to  the  changes  measured  by  conventional  MRI.  In  small  voxels,  the  detailed 
shape  of  the  water  resonance  is  likely  to  reflect  the  local  magnetic  environment  and 
changes  in  the  environment  caused  by  contrast  agents.  Each  resolvable  component  of  the 
water  resonance  may  represent  a  different  subvoxel  population  of  water  molecules  in  a 
discrete  environment.  Thus,  changes  in  resolvable  components  of  the  water  signal 
following  contrast  media  injection  may  reveal  details  of  subvoxel  physiology  and 
anatomy  that  would  be  difficult  or  impossible  to  detect  with  conventional  MRI.  For 
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example,  shoulders  that  appear  on  the  water  resonance  following  contrast  media  injection 
may  reflect  the  presence  of  subvoxelar  regions  with  dense  vasculature  -  and  this  may  be 
an  important  marker  for  cancer.  Thus  the  highest  spectral  resolution  possible,  given 
limitations  of  SNR  and  image  acquisition  time,  may  provide  an  important  view  of  the 
tumor  micro-environment.  Increases  in  signal-to-noise  ratio  due  to  improved  coil  design 
or  higher  field  would  allow  higher  spectral  resolution  in  the  future. 

Maximum  resolution  is  not  optimal  for  fat  suppression.  Rather,  the  FSR  is  maximized  at 
approximately  N  =  32,  after  which  there  is  a  slight  decrease  in  FSR.  Presumably,  this  is 
due  to  the  reduced  signal-to-noise  ratio  at  higher  spectral  resolutions.  Therefore,  for 
optimal  water  peak  height  images,  under  the  present  experimental  conditions,  it  may  be 
desirable  to  smooth  the  data  to  a  resolution  of  5  - 10  Hz,  while  optimal  detection  of 
effects  of  contrast  agents  requires  higher  resolution. 

Fitting  the  fat  resonance  to  a  Lorentzian  function  significantly  improves  fat  suppression 
(Fig.  2a),  allowing  detection  of  small  amounts  of  water  in  voxels  where  the  fat  signal  is 
predominant.  This  may  increase  sensitivity  to  small  and/or  diffuse  cancers.  The 
improvement  with  fitting  of  the  fat  resonance  comes  from  removal  of  both  the  wings  of 
signal  that  may  obscure  a  small  water  resonance,  and  the  baseline  (noise)  on  which  the 
water  signal  lies.  It  is  not  clear  that  equivalent  results  could  be  achieved  with  fat 
saturation.  If  the  fat  resonance  is  homogeneously  broadened  -  then  efficient  saturation  at 
the  center  of  the  resonance  would  reduce  the  amplitude  of  the  fat  signal  all  across  the 
spectrum.  However,  our  results  point  to  inhomogeneous  broadening:  Fig.  2b  shows  off- 
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resonance  fat  signal  of  magnitude  comparable  to  the  water  signal.  This  is  consistent  with 
NMR  studies  of  biopsy  samples,  [39-42]  demonstrating  that  the  fat  signal  represents 
multiple  chemical  species  with  varying  chemical  shifts.  In  that  case,  we  can  expect  high 
spectral  resolution  data  to  yield  superior  fat  suppression. 

In  conclusion,  we  have  demonstrated  that  at  1.5  Tesla  improvements  in  fat  suppression 
can  be  obtained  at  modest  spectral  resolution  (~10  Hz).  However,  improvements  in 
anatomic  detail  in  tumors,  and  increased  sensitivity  to  contrast  agents,  require  higher 
spectral  resolution.  This  higher  spectral  resolution,  combined  with  high  spatial  resolution 
may  improve  early  and  accurate  detection  of  breast  cancer. 
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Captions  for  Illustrations 

Figure  1:  FSR(a)and  Sr  (b),  calculated  at  varying  N  (varying  frequency  resolution), 
are  shown  for  five  patients  with  lesions  in  breast.  Crosses  represent  values  obtained 
without  fitting  the  fat  resonance,  while  the  diamonds  represent  the  corrected  values 
obtained  using  a  Lorentzian  fit  to  the  fat  resonance.  For  both  FSR  and  Sr,  the  fitting 
improves  the  quality  of  image. 

Figure  2:  a)  A  proton  spectrum  for  a  representative  voxel  containing  primarily  fat  (solid 
line)  and  the  Lorentzian  fit  for  the  fat  peak  (dashed  line)  are  shown,  b)  The  same 
spectrum  after  the  Lorentzian  fit  for  the  fat  signal  has  been  subtracted  (solid  line)  shows 
the  water  peak  much  more  clearly.  A  Lorentzian  fit  to  the  water  peak  is  shown  (dashed 
line).  Off-resonance  components  of  the  fat  peak  are  also  evident  at  ~  4  ppm.  The  vertical 
axis  is  shown  in  arbitrary  units. 

Figure  3:  The  water  (left  column)  and  fat  (right  column)  peak  height  images  of  a  patient 
with  infiltrating  DCIS,  calculated  using  varying  echo  train  lengths:  N  =  128  (a);  32  (b); 
and  8  (c)  are  shown.  There  is  improved  contrast  within  the  lesion  for  higher  N  and 
corresponding  spectral  resolution. 

Figure  4:  The  proton  spectra  in  five  sample  voxels  from  one  patient  with  infiltrating 
ductal  carcinoma  (pixels  #1-3)  and  one  healthy  volunteer  (pixels  #  4,  5),  measured 
using  varying  echo  train  lengths:  N  —  128  (a);  32  (b);  and  8(c)  are  shown.  The  arrow 
points  to  shoulders  in  the  water  resonance  which  are  not  discernible  in  lower  spectral 
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resolution  images.  The  water  peak  is  placed  at  0  ppm  in  all  spectra.  The  vertical  axis  is 
shown  in  arbitrary  units. 

Figure  5:  The  pre-  (solid  line)  and  post-contrast  (dashed  line)  spectra  from  a  single  voxel 
are  shown  at  the  maximum,  N  =  128  (a),  and  a  lower,  N  =  32  (b)  resolution.  The  arrow 
points  to  a  shoulder  on  the  water  resonance,  which  is  better  resolved  post-contrast  at  high 
resolution,  but  essentially  unchanged  in  low  resolution.  The  water  peak  is  placed  at  0 
ppm. 

Figure  6:  The  inhomogeneously  broadened  pixels  are  highlighted  for  images  obtained 
with  N  =  128  (a),  64  (b)  and  32  (c). 
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Breast  MR  Imaging  with 
High  Spectral  and  Spatial 
Resolutions:  Preliminary 
Experience1 


The  authors  evaluated  magnetic  reso¬ 
nance  (MR)  imaging  with  high  spec¬ 
tral  and  spatial  resolutions  (HSSR)  of 
water  and  fat  in  breasts  of  healthy  vol¬ 
unteers  (n  =  6)  and  women  with  sus¬ 
picious  lesions  (n  =  6).  Fat  suppression, 
edge  delineation,  and  image  texture 
were  improved  on  MR  images  derived 
from  HSSR  data  compared  with  those 
on  conventional  MR  images.  HSSR  MR 
imaging  data  acquired  before  and  af¬ 
ter  contrast  medium  injection  showed 
spectrally  inhomogeneous  changes 
in  the  water  resonances  in  small  vox¬ 
els  that  were  not  detectable  with 
conventional  MR  imaging. 

°  RSNA,  2002 


In  this  article,  we  describe  the  use  of 
echo-planar  spectroscopic  (EPS)  imaging 
for  magnetic  resonance  (MR)  imaging  of 
the  breast  with  high  spectral  and  spatial 
resolutions  (HSSR).  The  goal  of  this  ap¬ 
proach  is  to  acquire  data  with  the  spatial 
resolution  of  typical  anatomic  images 
and  with  a  high-resolution  spectrum  of 
the  water  and  fat  signals  for  each  image 
pixel.  Images  calculated  from  these  data 
sets  often  provide  improved  image  con¬ 
trast,  anatomic  accuracy,  and  functional 
information  (1—8).  HSSR  may  be  particu¬ 
larly  useful  for  MR  imaging  of  the  breast 
because  they  may  improve  separation  of 
water  and  fat  signals,  contrast,  edge  de¬ 
lineation,  and  sensitivity  to  contrast  me¬ 
dia  and,  thus,  increase  the  sensitivity  for 
early  carcinomas. 

The  present  application  of  HSSR  MR 
imaging  is  a  natural  extension  of  previ¬ 
ous  work  that  demonstrated  the  impor¬ 
tance  of  incorporating  spectroscopic  in¬ 
formation  into  MR  images.  Dixon  first 
demonstrated  that  two  points  of  spectral 
resolution  often  allow  good  separation  of 


water  and  fat  signals  (9).  Glover  and 
Schneider  (10),  Glover  (11),  and  others 
(12,13)  produced  further  improvements 
in  fat-water  separation  and  correction  for 
B0  inhomogeneity  with  several  points  of 
spectral  resolution.  More  recently,  ad¬ 
vances  in  MR  hardware  and  software 
have  allowed  rapid  acquisition  of  data 
with  use  of  EPS  MR  imaging  (14,15)  at 
very  high  spectral  and  spatial  resolutions. 

With  EPS  MR  imaging,  details  of  the 
water  and  fat  line  shapes  in  each  small 
image  pixel  can  be  resolved  (4,16)  with 
reasonable  acquisition  times.  HSSR  MR 
imaging  is  important  when  the  water  and 
fat  resonances  in  small  pixels  have  com¬ 
plicated  line  shapes.  Inhomogeneous 
broadening  due  to  local  magnetic  suscep¬ 
tibility  gradients  is  expected  on  theoretic 
grounds  (17-20)  and  has  been  demon¬ 
strated  experimentally  at  high  spatial  res¬ 
olution  (1,3,4,21-23).  This  effect  is  espe¬ 
cially  strong  in  and  near  tumors,  where 
deoxygenated  blood  and  a  generally  het¬ 
erogeneous  environment  produce  broad 
complicated  water  lines  (1,22).  The  de¬ 
tailed  structure  of  the  water  line  can  be 
an  important  source  of  information  re¬ 
garding  local  (subpixel)  anatomy  and 
physiology  and  can  be  analyzed  to  im¬ 
prove  contrast  and  accuracy  on  HSSR  MR 
images  (1,24). 

The  EPS  MR  imaging  approach  origi¬ 
nally  developed  by  Doyle  and  Mansfield 
(14)  and  Mansfield  (15)  and  subsequently 
refined  by  others  (25)  allows  rapid  acqui¬ 
sition  of  HSSR  MR  imaging  data  sets.  For 
example,  256  x  256  matrices  in  the  spa¬ 
tial  domain  combined  with  128-256 
points  in  the  spectral  dimension  and 
2-Hz  spectral  resolution  can  be  acquired 
in  less  than  2  minutes  with  the  simplest 
k-space  sampling  methods.  Much  more 
rapid  data  acquisition  is  possible  when  k 
space  is  sampled  along  two  or  three  axes 
during  readout.  Results  to  date  suggest 
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that  the  rapidly  switching  gradients  re¬ 
quired  by  EPS  MR  imaging  on  clinical  MR 
imagers  produce  relatively  little  distor¬ 
tion  in  proton  line  shapes,  even  at  reso¬ 
lutions  of  1-2  Hz.  These  results  suggest 
that  the  effects  of  eddy  currents  are  rela¬ 
tively  small  when  state-of-the-art  self- 
shielded  gradient  coils  are  used  (4,6). 
Thus,  it  is  realistic  to  consider  application 
of  HSSR  MR  imaging  to  evaluate  suspi¬ 
cious  breast  lesions  in  a  clinical  setting. 
The  purposes  of  this  study  were  to  per¬ 
form  a  preliminary  evaluation  of  HSSR 
MR  imaging  of  the  breast  and  suspicious 
breast  lesions  and  to  compare  HSSR  MR 
images  with  conventional  T2-weighted 
and  T1 -weighted  fat-saturated  MR  im¬ 
ages. 


I  Materials  and  Methods 

Data  Acquisition 

Six  healthy  female  volunteers  (age 
range,  28-43  years;  mean  age,  37  years  ± 
5  [mean  ±  SD])  and  six  female  patients 
(age  range,  40-61  years;  mean  age,  47 
years  ±  7)  were  enrolled  in  this  study. 
The  healthy  volunteers  (ie,  women  with 
no  history  of  abnormalities  at  mammog¬ 
raphy  or  palpation)  were  recruited  from 
among  the  staff  of  the  radiology  depart¬ 
ment.  The  patients  were  recruited 
through  the  breast  clinic  on  the  basis  of 
having  suspicious  lesions  with  a  largest 
diameter  of  4  mm  to  3  cm  on  mammo¬ 
grams.  All  six  patients  were  considered  to 
have  high  risk  of  breast  cancer.  Three 
patients  (selected  by  U.B.)  underwent 
breast  biopsy  (performed  by  the  clinical 
pathology  service);  one  of  them  was 
found  to  have  a  benign  lesion,  and  the 
other  two  were  determined  to  have  infil¬ 
trating  ductal  carcinomas.  The  other 
three  patients  did  not  undergo  biopsy.  In 
three  of  the  six  patients,  abnormalities 
were  not  depicted  clearly  on  the  conven¬ 
tional  MR  images,  probably  because  only 
a  limited  number  of  sections  were  im¬ 
aged.  MR  imaging  in  all  subjects  was  per¬ 
formed  according  to  a  protocol  approved 
by  the  institutional  review  board  after 
informed  consent  had  been  obtained. 

The  patients  underwent  MR  imaging 
before  and  after  intravenous  injection  of 
a  gadolinium  chelate  (Omniscan;  Ny- 
comed  Amersham,  Princeton,  NJ)  at  a 
dose  of  0.2  mmol  per  kilogram  of  body 
weight.  Approximately  20  mL  of  the  con¬ 
trast  media  solution  was  injected  by  us¬ 
ing  an  automatic  injector  (Medrad,  Pitts¬ 
burgh,  Pa)  at  a  rate  of  2  mL/sec.  The 
healthy  volunteers  did  not  receive  con¬ 
trast  media. 


The  MR  images  were  acquired  with  a 
1.5-T  clinical  system  (Signa;  GE  Medical 
Systems,  Milwaukee,  Wis)  equipped  with 
gradients  with  a  maximum  slew  rate  of 
120  mT/m/sec  and  maximum  amplitude 
of  23  mT/m.  Signal  was  detected  by  using 
a  dedicated  phased-array  breast  coil,  and 
shimming  was  performed  by  using  the 
signal  from  the  entire  sensitive  volume  of 
the  coil.  HSSR  MR  images  were  acquired 
with  EPS  MR  imaging  sequences.  After 
section  selection  (section  thickness,  4 
mm)  and  phase  encoding  (256  steps), 
128  gradient  echoes  were  acquired  by  us¬ 
ing  trapezoidal  gradient  pulses  with  alter¬ 
nating  polarity.  A  crusher  gradient  was 
applied  at  the  end  of  the  echo  train  to 
eliminate  artifacts  due  to  residual  trans¬ 
verse  magnetization.  Each  gradient  echo 
was  sampled  at  256  points.  The  data  were 
digitized  at  a  bandwidth  of  plus  or  minus 
62.5  kHz,  and  the  time  between  the  cen¬ 
ters  of  gradient  echoes  was  3.0  msec.  The 
proton  free  induction  decay  was  sampled 
for  384  msec,  and  the  repetition  time  was 
500  msec.  The  resulting  data  had  an  in¬ 
plane  spatial  resolution  of  less  than  1 
mm  (field  of  view,  <24  cm),  spectral  res¬ 
olution  of  2.6  Hz,  and  spectral  band¬ 
width  of  333  Hz.  This  bandwidth  was 
sufficient  to  resolve  the  water  and  fat  res¬ 
onances  at  1.5  T.  Some  spectral  wrap¬ 
around  was  corrected.  Both  spectral  and 
spatial  resolutions  were  sufficient  to 
avoid  truncation  artifacts.  Sagittal  sec¬ 
tions  were  imaged  with  the  readout  gra¬ 
dient  applied  in  the  anteroposterior  di¬ 
rection  to  minimize  artifacts  due  to 
respiratory  and  cardiac  motion. 

A  T2-weighted  multisection  fast  spin- 
echo  sequence  was  used  to  determine  lo¬ 
cations  of  suspicious  lesions  (repetition 
time  msec/echo  time  msec  [effective]  of 
4,000/176;  echo  train  length,  12;  section 
thickness,  4  mm;  matrix,  256  X  256;  field 
of  view,  ~24  mm;  two  signals  acquired). 
T1 -weighted  MR  images  with  fat  satura¬ 
tion  were  acquired  with  either  the  fast 
spin-echo  method  (600/14  [effective]; 
echo  train  length,  12;  section  thickness,  4 
mm;  matrix,  256  X  256;  field  of  view, 
~24  mm;  eight  signals  acquired)  or  the 
conventional  spin-echo  method  (617/14; 
section  thickness,  4  mm;  matrix,  256  x 
256;  field  of  view,  —24  cm;  one  signal 
acquired).  HSSR  MR  images  of  the  se¬ 
lected  section  were  acquired  as  described 
earlier,  that  is,  two  images  before  contrast 
media  injection  and  one  image  between 
1  and  3  minutes  after  injection.  The  ac¬ 
quisition  time  for  one  HSSR  MR  image 
was  128  seconds.  Tl-weighted  fat-satu¬ 
rated  MR  imaging  was  repeated  at  3  min¬ 
utes  after  contrast  media  injection. 


To  evaluate  the  performance  of  the  EPS 
MR  imaging  pulse  sequence,  spherical 
(-10-cm-diameter)  water  phantoms  con¬ 
taining  10  mmol/L  copper  sulfate  were 
imaged.  The  water  signal  in  small  voxels 
in  these  phantoms  (as  measured  with  sin¬ 
gle-voxel  spectroscopy)  had  a  line  width 
of  approximately  1-2  Hz.  The  line  width 
in  small  voxels  of  high-resolution  spec¬ 
troscopic  images  of  the  phantoms  ac¬ 
quired  with  EPS  MR  imaging  was  an  ac¬ 
curate  indicator  of  the  eddy  currents 
produced  by  the  EPS  MR  imaging  pulse 
sequence.  Eddy  currents  would  signifi¬ 
cantly  broaden  the  lines  and  cause  phas¬ 
ing  artifacts  due  to  frequency  shifts  dur¬ 
ing  the  decay  of  the  free  induction  decay. 
EPS  MR  images  were  acquired  with  the 
same  parameters  as  were  used  for  breast 
imaging,  as  described  earlier. 

Data  Analysis  and  Synthesis  of  Images 

Data  analysis  was  performed  with  in- 
house  software  written  in  IDL  (Interactive 
Data  Language;  Research  Systems,  Boulder, 
Colo).  Raw  data  generated  with  EPS  MR 
imaging  were  processed  by  using  methods 
described  elsewhere  (4).  In  summary,  a 
three-dimensional  Fourier  transform  with 
respect  to  two  k-space  axes  and  the  evolu¬ 
tion  of  the  free  induction  decay  provide  a 
data  set  in  which  spectral  and  spatial  infor¬ 
mation  are  well  separated.  A  high-resolu¬ 
tion  proton  spectrum  is  associated  with 
each  pixel  in  the  image.  After  the  water 
and  fat  signals  are  separated,  a  variety  of 
images  based  on  the  water  and  fat  line 
shapes  can  be  produced  (1,4,8,22,26). 

For  the  purposes  of  the  present  study, 
images  were  synthesized  with  signal  in¬ 
tensity  proportional  to  the  signal  peak 
heights  of  water  and  fat  (for  mixed  T2* 
and  T1  contrast),  integrals  (for  T1  con¬ 
trast),  line  widths  (for  T2*  contrast),  and 
resonance  frequencies.  In  addition,  im¬ 
ages  were  produced  to  highlight  pixels 
with  highly  asymmetric  water  peaks  or 
with  multiple  resolved  components  of 
the  water  resonance.  To  represent  the  ef¬ 
fects  of  contrast  media,  images  were  cal¬ 
culated  to  show  changes  in  peak  height, 
frequency  shifts,  and  spectrally  heteroge¬ 
neous  effects. 

Separation  of  Water  and  Fat  Signals 

Nonuniformity  of  the  magnetic  field 
across  the  breast  caused  variations  in  the 
resonance  frequencies  of  water  and  fat. 
To  identify  water  and  fat  signals  in  each 
pixel,  a  frequency  map  was  constructed 
by  using  the  following  algorithm. 

1.  The  largest  spectral  peak  was  auto¬ 
matically  identified  in  all  image  pixels 


578  *  Radiology  •  August  2002 


Du  et  al 


with  signal  intensities  three  times  greater 
than  the  noise  level.  A  rough  frequency 
map  was  constructed  to  represent  the  res¬ 
onance  frequencies  of  the  identified 
peaks.  This  map  showed  regions  that  pri¬ 
marily  contained  water  (hereafter,  a  wa¬ 
ter  region)  and  regions  that  primarily 
contained  fat  (hereafter,  a  fat  region). 
The  water  and  fat  regions  composed  most 
but  not  all  of  the  breast.  Some  isolated 
pixels  were  not  included  in  either  region. 
Abrupt  transitions  in  resonance  fre¬ 
quency  of  the  largest  spectral  peak  (—220 
Hz)  occurred  at  the  interfaces  between 
water  and  fat  regions.  Other  sharp  tran¬ 
sitions  equal  to  the  spectral  bandwidth 
(333  Hz)  occurred  as  a  result  of  fold  back. 
Within  each  region,  the  frequencies  of 
water  or  fat  were  assumed  to  vary  slowly 
(ie,  there  were  no  sharp  magnetic  field 
gradients). 

2.  A  seed  pixel  was  manually  selected 
(W.D,  and  G.S.K.  together)  in  each  of  one 
to  three  fat  regions  and  in  each  of  one  to 
three  water  regions.  Water  or  fat  signals 
were  identified  in  pixels  adjacent  to  the 
seed  pixel  by  searching  within  a  narrow 
bandwidth  around  the  water  or  fat  reso¬ 
nance  frequencies  in  the  seed  pixel.  This 
process  was  implemented  iteratively  by 
using  a  region-growing  program  written 
(W.D.)  in  IDL.  Within  each  fat  region, 
the  water  resonance  frequency  was  calcu¬ 
lated  in  each  pixel  by  using  the  known 
fat  resonance  frequency  and  the  chemi¬ 
cal  shift  offset  (+220  Hz).  Similarly,  the 
fat  resonance  frequency  was  determined 
in  water  regions.  Fold-back  effects  were 
identified  on  the  basis  of  the  sharp 
change  in  water  or  fat  resonance  frequen¬ 
cies  and  were  corrected  by  adding  or  sub¬ 
tracting  the  spectral  bandwidth  (333  Hz, 
or  128  frequency  bins)  to  the  frequency 
of  the  folded  resonance.  As  a  result  of  this 
correction,  the  resonance  frequencies  for 
water  and  fat  varied  smoothly  across  all 
fat  and  water  regions. 

3,  In  pixels  outside  the  fat  and  water 
regions — usually  pixels  with  low  signal- 
to-noise  ratios — the  water  and  fat  reso¬ 
nance  frequencies  were  estimated  from 
the  frequency  values  in  the  water  or  fat 
regions  by  using  a  morphologic  dilation 
operation.  The  gray-value  dilation  opera¬ 
tion  was  repeatedly  applied  on  the  fre¬ 
quency  map  until  the  gaps  between  the 
water  or  fat  regions  were  filled,  and  the 
water  and  fat  resonance  frequencies  in  all 
pixels  throughout  the  breast  were  deter¬ 
mined. 

The  water  spectrum  was  taken  to  be 
the  signal  within  plus  or  minus  10  fre¬ 
quency  bins  around  the  central  water  fre¬ 
quency,  as  obtained  from  the  frequency 


map.  To  calculate  water  signal  peak 
height  and  integral,  the  water  line  was  fit 
to  the  magnitude  of  a  Lorentzian  line 
plus  a  linear  baseline  by  using  a  nonlin¬ 
ear  least  squares  curve-fitting  algorithm. 
The  water  peak  height  was  taken  as  the 
highest  point  of  the  water  line  after  the 
baseline  was  removed  with  the  spectral 
fitting  routine.  In  principle,  this  proce¬ 
dure  could  not  precisely  fit  water  lines  in 
many  pixels,  particularly  those  that  were 
strongly  asymmetric  or  contained  multiple 
resolved  components  (23).  In  practice, 
however,  the  fitting  procedure  provided 
good  approximations  for  calculation  of 
peak-height  and  integral  images.  The 
procedure  allowed  accurate  detection  of 
a  small  amount  of  water,  even  In  the 
presence  of  a  sloping  baseline  due,  for 
example,  to  large  fat  signals  with  broad 
wings  or  signals  that  leaked  from  other 
voxels. 

Texture  Analysis 

Image  texture  was  measured  In  a  re¬ 
gion  of  interest  (ROI)  selected  to  include 
as  much  as  possible  of  the  breast  but  to 
exclude  the  chest  wall.  ROIs  were  se¬ 
lected  manually  by  MR  physicists  (X.F. 
and  W.D.  together)  with  advice  from  a 
radiologist  (P.M.M.  or  U.B.).  Image  signal 
intensity  In  the  ROI  was  normalized  by 
the  average  value  of  image  signal  inten¬ 
sity  in  the  ROI  so  that  images  from  dif¬ 
ferent  examinations  could  be  compared. 
Image  signal  intensity  was  expressed  as  a 
function  (f)  of  the  column  (x)  and  row  (y) 
index  for  each  pixel.  Mathematically,  f(x, 
y)  represents  a  surface  spanning  an  ROI 
on  an  image  coordinate  plane  x-y.  The 
area  of  the  surface,  5,  was  found  from 

S  =  S^T+fI+ff  (1) 
where  the  summation  is  over  the  ROI, 
4  =  8/78*,  and  4  =  8 ffhy.  We  defined  the 
normalized  surface  area  (Sr)  as 

Sr  =  100.0  •  I S(ROI)  -  Sxy  I/S  ' y  (2) 

where  S(ROI)  is  the  surface  area  calcu¬ 
lated  from  Equation  (1)  and  Sxy  is  the  area 
of  S(ROI)  projected  onto  the  x-y  plane. 
Thus,  Sr  is  a  measure  of  the  variation  or 
roughness  in  image  signal  intensity  (eg,  if 
the  image  signal  intensity  is  constant, 
then  Sr  is  0.0). 

Locating  Pixels  with  Multiple  Peaks 
and  Asymmetric  Peaks 

Multiple  peaks  and  asymmetric  peaks 
were  identified  automatically  with  soft¬ 
ware  written  in  IDL.  To  find  the  number 
of  resolved  peaks  in  the  water  resonance 


in  each  pixel,  all  possible  peaks  in  the 
water  line  were  automatically  located 
with  the  software  (ie,  those  points  in  the 
spectrum  where  the  signal  was  larger 
than  that  at  the  adjacent  points).  Then, 
with  the  program,  spectral  amplitude  was 
measured  at  the  local  minima  between 
these  candidate  peaks.  Candidate  peaks 
were  selected  as  true  peaks  if  the  differ¬ 
ence  between  the  local  maximum  and 
minimum  was  two  times  larger  than  the 
root-mean-square  of  the  noise  level.  The 
number  of  resolved  peaks  was  counted 
for  each  pixel,  and  images  were  con¬ 
structed  to  show  the  location  of  pixels 
that  contained  multiple  water  peaks. 

To  determine  the  asymmetry  of  mag¬ 
nitude  spectra,  the  sides  of  the  spectra  on 
either  side  of  the  global  maxima  were 
compared.  Spectra  were  identified  as 
asymmetric  if  the  difference  between  the 
integrals  of  the  high-  and  low-frequency 
sides  of  the  resonance  was  larger  than 
10%  of  the  total  integral. 

Analysis  of  Effects  of  Contrast  Media 
on  HSSR  Data 

Before  the  effects  of  contrast  media 
were  evaluated,  the  HSSR  data  sets  were 
corrected  for  motion  in  two  dimensions. 
In-plane  breast  motions,  including  trans¬ 
lation,  rotation,  and  scaling,  were  mod¬ 
eled  with  a  six-parameter  coordinate 
transform.  The  postcontrast  peak-height 
image  was  taken  as  the  reference  to 
which  the  precontrast  peak-height  image 
was  registered.  Six  to  12  corresponding 
points  that  represented  shared  anatomic 
features  were  manually  identified  on  the 
two  peak-height  images.  The  parameters 
of  the  coordinate  transform  were  deter¬ 
mined  by  minimizing  the  squared  differ¬ 
ence  between  the  coordinates  of  the  cor¬ 
responding  points.  Finally  the  coordinate 
transform  and  bilinear  interpolation 
were  applied  to  obtain  the  motion-cor¬ 
rected  precontrast  HSSR  data  set. 

Software  was  written  in  IDL  to  measure 
the  change  in  the  water  signal  integral, 
peak  height,  frequency,  and  T2*  caused 
by  the  contrast  media.  The  integral  of  the 
water  magnitude  spectrum  before  con¬ 
trast  media  injection  was  subtracted  from 
that  after  contrast  media  injection.  The 
frequency  shift  due  to  the  contrast  media 
in  each  pixel  was  calculated  after  zero 
filling  of  the  free  induction  decays  to 
2,048  points.  T2*  was  estimated  by  divid¬ 
ing  the  peak  height  by  the  peak  integral 
of  the  water  spectrum  for  each  pixel.  The 
changes  in  T2*  due  to  the  contrast  media 
were  then  computed. 

To  study  the  spectrally  inhomogeneous 
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effects  of  the  contrast  media  (1,3, 4,8, 
22,26),  images  were  generated  on  the  ba¬ 
sis  of  the  difference  spectrum  obtained 
by  subtracting  the  precontrast  spectrum 
in  each  pixel  from  the  postcontrast  spec¬ 
trum.  Both  the  pre-  and  postcontrast  free 
induction  decays  were  zero  filled  to  512 
points  before  being  Fourier  transformed. 
The  postcontrast  spectrum  from  each 
pixel  was  then  shifted  until  the  absolute 
difference  between  the  pre-  and  postcon¬ 
trast  spectra  was  minimized.  This  shift  cor¬ 
rected  for  contrast  media-induced  changes 
in  bulk  susceptibility  that  shifted  the  en¬ 
tire  water  signal.  The  resulting  difference 
spectra  were  used  to  synthesize  images 
on  which  signal  intensity  was  propor¬ 
tional  to  the  greatest  positive  and  nega¬ 
tive  amplitudes  in  the  difference  spectra 
and  also  to  the  greatest  absolute  change 
in  spectral  amplitude  in  each  pixel.  In 
addition,  we  synthesized  images  propor¬ 
tional  to  the  product  of  absolute  values 
of  the  greatest  positive  and  negative 
changes  in  the  difference  spectrum. 

Evaluation  by  Experienced 
Radiologists 

The  conventional  and  HSSR  MR  im¬ 
ages  of  the  breast  in  all  subjects  (n  =  12) 
were  reviewed  independently  by  two  ra¬ 
diologists  (U.B.  and  who  have  10 

and  5  years  of  experience,  respectively). 
The  radiologists  were  provided  with  the 
medical  histories  including  the  mammo- 
graphic  findings,  if  any,  but  excluding 
the  biopsy  results.  The  radiologists  eval¬ 
uated  the  diagnostic  usefulness  of  the 
HSSR  MR  images  over  that  of  the  conven¬ 
tional  images  in  terms  of  the  separation 
of  fat  and  water  signals,  the  depiction  of 
normal  tissues,  lesions,  and  surrounding 
vascular  structures,  and  the  sensitivity  of 
the  imaging  technique  to  the  contrast 
media. 

Statistical  Analysis 

Nonparametric  (Wilcoxon  signed- 
rank)  tests  were  performed  with  commer¬ 
cially  available  software  (Minitab,  State 
College,  Pa)  with  data  from  the  normal 
breasts.  The  tests  were  performed  to  com¬ 
pare  the  effectiveness  of  fat  suppression 
and  image  texture  or  sharpness  on  water 
HSSR  MR  images  and  conventional  Tl- 
weighted  fat-saturated  MR  images.  Differ¬ 
ences  with  P  <  .05  were  considered  sta¬ 
tistically  significant. 

I  Results 

Accurate  high-resolution  spectroscopic 
imaging  with  EPS  MR  imaging  requires 


that  the  gradient  waveforms  be  produced 
with  minimal  eddy  currents  and  that 
shimming  is  adequate.  EPS  MR  imaging 
of  water  phantoms  produced  narrow 
(<2-Hz)  symmetric  water  lines  with  no 
phasing  artifacts.  These  images  demon¬ 
strated  that  eddy  current  and  shimming 
effects  were  small  and  did  not  distort 
spectra,  even  at  resolution  of  less  than 
2  Hz. 

The  high-resolution  EPS  MR  images 
demonstrated  strong  spectral  inhomoge¬ 
neity  of  water  signal  line  shapes  in  the 
breast.  Water  signals  from  individual 
small  voxels  often  contained  multiple  re¬ 
solved  components  or  were  clearly  asym¬ 
metric.  For  example,  the  water  resonance 
in  12.5%  ±  9.0  (mean  ±  SD)  of  all  pixels 
in  six  normal  breasts  was  composed  of 
two  or  more  multiple  resolved  compo¬ 
nents.  A  larger  percentage  of  pixels  in 
these  breasts  contained  asymmetric  pix¬ 
els:  20.4%  ±  10.8.  These  parameters  were 
similar  in  the  group  of  women  with 
breast  lesions:  20.7%  ±  8.2  of  pixels  con¬ 
tained  water  resonances  with  multiple  re¬ 
solved  components  and  30.7%  ±  10.9 
contained  asymmetric  resonances. 

These  percentages  are  based  on  the  to¬ 
tal  number  of  pixels  in  the  breast  and, 
thus,  are  low  owing  to  the  fact  that  some 
pixels  do  not  contain  detectable  water. 
In  pixels  that  contained  a  significant 
amount  of  MR-detectable  water,  approx¬ 
imately  30.4%  ±  6.3  of  the  water  reso¬ 
nances  contained  multiple  resolved  com¬ 
ponents.  After  contrast  media  injection 
(n  -  6),  16.9%  ±  8.8  of  pixels  in  the 
breasts  contained  water  signals  with  mul¬ 
tiple  resolved  components.  The  spectral 
inhomogeneity  of  the  water  resonance  is 
illustrated  in  Figure  1.  An  image  synthe¬ 
sized  from  water  signal  peak  height  and 
representative  spectra  acquired  before 
and  after  contrast  media  injection  from 
selected  pixels  is  shown.  Inhomogeneous 
broadening  of  the  water  peak  is  clear  in 
many  of  the  spectra,  and  some  water  res¬ 
onances  have  two  or  more  resolvable 
components  (spectra  3,  4,  5,  and  7).  The 
response  of  the  water  resonance  to  the 
contrast  media  is  different  from  that  of 
the  fat  resonance.  In  spectrum  1,  for  ex¬ 
ample,  the  water  peak  shifts  to  higher 
frequency  and  its  magnitude  increases, 
while  the  fat  peak  does  not  shift  and  its 
magnitude  decreases  slightly.  In  spectra 
4,  5,  and  7,  the  contrast  agent  produces  a 
well-defined  shoulder  or  fully  resolved 
component  on  the  high-frequency  side 
(low-field-strength  side)  of  water.  The  in¬ 
creases  in  overall  peak  height  of  the  water 
signal  are  expected  because  the  acquisi¬ 
tions  are  T1  weighted. 


Inhomogeneous  broadening  of  the  wa¬ 
ter  resonance  can  be  used  to  produce  im¬ 
ages  that  highlight  anatomic  features. 
For  one  of  the  breasts  imaged,  the  spatial 
distribution  of  pixels  with  non-Lorent- 
zian  water  spectra  is  illustrated  in  Figure 
2,  which  highlights  pixels  in  which  the 
water  resonance  has  multiple  compo¬ 
nents  (Fig  2,  B)  or  is  asymmetric  (Fig  2, 
C).  In  this  breast,  the  water  resonances  in 
approximately  21%  of  the  tumor  pixels 
contain  two  or  more  resolved  compo¬ 
nents,  and  the  spectra  in  39%  of  pixels 
are  asymmetric.  Pixels  with  complicated 
water  resonances  tended  to  be  near  the 
tumor  edges  or  near  large  vessels. 

Although  asymmetry  and  multiple  re¬ 
solved  components  are  clear  indicators  of 
inhomogeneous  broadening,  even  rela¬ 
tively  symmetric  water  resonances  are  of¬ 
ten  inhomogeneously  broadened.  Spec¬ 
tral  inhomogeneity  is  clear  from  the 
effects  of  contrast  media  on  distinct  com¬ 
ponents  of  relatively  symmetric  water 
resonances  (eg,  spectra  3,  4,  5,  and  7  in 
Fig  1).  The  effects  of  contrast  media  in 
many  pixels  were  spectrally  inhomoge¬ 
neous,  that  is,  the  changes  caused  by  the 
contrast  media  were  not  uniform  across 
the  water  spectrum.  The  largest  change 
in  the  spectrum  was  often  not  at  the  cen¬ 
ter  of  the  water  resonance. 

To  evaluate  fat  suppression  on  water 
HSSR  MR  images,  the  signal  in  parenchy¬ 
mal  pixels  on  these  images  (ie,  pixels  that 
primarily  contained  water)  was  divided 
by  the  signal  in  pixels  that  primarily  con¬ 
tained  fat.  This  ratio  is  referred  to  as  the 
fat-suppression  ratio  in  the  following 
data.  The  fat-suppression  ratios  were  sig¬ 
nificantly  higher  (Wilcoxon  signed-rank 
test,  two-tail  P  =  .036)  on  water  HSSR  MR 
images  (3.7  ±  0.1  [  mean  ±  standard  er¬ 
ror  of  the  mean])  than  on  Tl-weighted 
fat-saturated  images  (1.4  ±  0.2).  Regions 
of  insufficient  fat  suppression  on  the 
conventional  images  were  evident,  but  in 
addition  the  general  background  in  pre¬ 
dominantly  fat  regions  on  conventional 
fat-saturated  MR  images  was  brighter 
than  that  on  water  HSSR  MR  images.  This 
difference  is  illustrated  in  Figure  3,  which 
shows  conventional  and  HSSR  MR  im¬ 
ages  of  single  sagittal  sections  in  six 
healthy  breasts. 

Improved  fat  suppression  on  water 
HSSR  MR  images  enhances  contrast  in 
the  parenchyma.  The  Table  shows  quan¬ 
titative  comparisons  of  the  Tl-weighted 
fat-saturated  MR  images  and  water  HSSR 
MR  images  presented  in  Figure  3.  Because 
the  signal-to-noise  ratios  were  similar  on 
HSSR  and  Tl-weighted  MR  images  (Ta¬ 
ble),  the  normalized  surface  area  (Sr  in  Eq 
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Figure  1*  Representative  spectra  from  selected  pixels  before  (dotted  lines)  and  after  (solid  lines)  contrast  medium  injection.  Reference 
image  shows  a  sagittal  section  through  a  breast  with  ductal  carcinomas  (sampled  with  spectra  1,  4,  5,  and  7).  All  spectra  are  on  the  same 
scale  and  are  referenced  to  the  same  carrier  frequency.  Frequency  shifts  due  to  B0  variation  across  the  breast  are  apparent.  The  water  peak 
(w)  and  fat  peak  (f)  are  labeled  for  each  spectrum.  Spectra  in  some  pixels  have  approximately  the  same  shape  before  and  after  contrast 
media  injection  (spectra  1,  2,  3,  and  6),  which  indicates  that  effects  of  interimage  motion  have  been  minimized. 


A  B  c 


Figure  2,  A,  Pixels  with  non-Lorentzian  spectral  shapes  are  overlaid  on  a  sagittal  reference  water 
peak-height  HSSR  MR  image  of  a  breast  with  a  lesion  (arrow).  B,  Pixels  in  which  the  water 
resonance  contains  multiple  resolved  peaks  are  highlighted.  Only  peaks  with  signal-to-noise  ratio , 
of  at  least  2  were  identified  as  secondary  peaks.  C,  Pixels  with  strongly  asymmetric  water 
resonances  are  highlighted. 


[2])  can  be  used  to  compare  texture  and 
sharpness  (Sr)  of  the  HSSR  MR  images 
with  those  of  the  T1 -weighted  MR  im¬ 
ages.  Sr  values,  tabulated  in  the  Table,  are 


larger  (Wilcoxon  signed-rank  test,  two- 
tailed  P  =  .036)  on  water  HSSR  MR  im¬ 
ages  than  those  on  Tl-weighted  MR  im¬ 
ages.  The  results  demonstrate  that  HSSR 


MR  images  generally  provide  better  fat- 
water  separation  and  better  contrast  and 
texture  compared  with  those  on  conven¬ 
tional  Tl-weighted  MR  images  with  fat 
saturation. 

Strong  enhancement  was  found  on  MR 
images  obtained  in  the  two  patients  with 
histopathologieally  proved  invasive  duc¬ 
tal  carcinomas.  Two  experienced  radiolo¬ 
gists  (blinded  independent  evaluation  by 
U.B.  and  P.M.M.)  agreed  that  image  con¬ 
trast  was  improved  on  HSSR  MR  images 
compared  with  that  on  conventional  im¬ 
ages  and  that  tumor  boundaries  were 
more  clearly  defined.  Fat  signal  was  more 
effectively  suppressed  on  the  water  HSSR 
MR  images  so  that  the  water-containing 
structures  were  clearly  seen  against  a  dark 
background.  These  findings  are  illus¬ 
trated  in  Figure  4,  where  conventional 
and  HSSR  MR  images  are  compared.  The 
contrast  ratio  for  the  first  lesion  (arrow7  in 
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part  A)  compared  with  surrounding  tis¬ 
sue  (average  signal  in  lesion  divided  by 
average  signal  in  surrounding  tissue)  is 
4.4  on  the  water  peak-height  HSSR  MR 
image  before  contrast  media  injection 
(part  C),  4.8  after  contrast  media  injec¬ 
tion  (part  D),  and  1.9  on  the  Tl-weighted 
fat-saturated  MR  images  after  injection 
(part  B).  The  vascular  network  (arrow  in 
part  D)  in  the  first  breast  is  better  de¬ 
picted  on  HSSR  MR  images.  Further,  U.B. 
and  P.M.M.  found  that  the  edges  of  the 
second  lesion  (arrow  in  part  E)  are  more 
clearly  delineated  on  peak-height  HSSR 
MR  images  than  on  conventional  MR  im¬ 
ages,  especially  near  the  chest  wall  (arrow 
in  part  H).  In  this  case,  HSSR  MR  images 
provide  a  clearer  view  of  whether  or  not 
the  lesion  is  infiltrating  the  chest  wall. 

The  effects  of  contrast  agents  on  water 
signal  line  shape  and  resonance  frequency 
can  be  used  to  produce  novel  images.  For 
example,  Figure  5  illustrates  the  effects  of 
contrast  media  on  HSSR  data  sets.  Figure  5, 
A  shows  changes  on  an  image  derived  from 
the  first  gradient  echo  of  the  EPS  echo  train 
(postcontrast  image  minus  precontrast  im¬ 
age).  This  image  simulates  the  effect  of 
the  contrast  media  on  conventional  Tl- 
weighted  gradient-echo  MR  images.  Strong 
signal  enhancement  is  seen  in  blood  ves¬ 
sels  at  the  tumor  rim  and  in  some  pixels 
inside  the  tumor.  Figure  5,  B  shows  the 
change  in  the  water  peak-height  image. 
The  inherent  T2*  weighting  in  this  image 
shows  regions  of  decreased  or  increased 
signal  intensity.  Figure  5,  C  shows  a  map  of 
the  water  peak-frequency  change  after  con¬ 
trast  media  injection.  These  changes  were 
greatest  in  large  vessels.  Frequency  changes 
are  evident  in  elongated  structures  in  the 
tumor  interior  that  perhaps  indicate  the 
primary  veins  and  arteries  in  the  tumor 
center.  Finally,  image  signal  intensity  in 
Figure  5,  D  is  proportional  to  the  product 
of  the  absolute  values  of  the  largest  positive 
and  negative  changes  in  the  water  spec¬ 
trum.  This  product  image  highlights  pixels 
with  spectrally  heterogeneous  contrast  me¬ 
dia  effects  and  shows  changes  beyond  the 
edge  of  the  upper  quadrant  of  the  lesion 
that  are  not  clear  on  the  conventional  im¬ 
age  (Fig  5,  A).  The  images  derived  from 
HSSR  data  show  strong  features  that  are 
not  evident  on  conventional  images  (Fig  4, 
A,  B). 


Figure  3.  Conventional  and  HSSR  MR  images  of  six  healthy  breasts.  For  each  breast  (one  row), 
four  sagittal  images  are  shown.  From  left  to  right:  T2-weighted  fast  spin-echo  MR  images, 
Tl-weighted  fat-saturated  fast  spin-echo  MR  images,  water  peak-height  HSSR  MR  images,  and  fat 
peak-height  HSSR  MR  images.  Image  signal  intensities  are  adjusted  to  a  similar  level,  and  images 
are  displayed  with  identical  window  settings.  Conventional  and  HSSR  MR  images  were  each 
acquired  in  approximately  1. 5-2.0  minutes. 


I  Discussion 

These  preliminary  results  demonstrate 
the  application  of  HSSR  imaging  to  hu¬ 
man  breast  MR  imaging  and  evaluation 
of  suspicious  lesions.  We  evaluated  both 


conventional  properties  (ie,  contrast,  tex¬ 
ture,  and  fat  suppression)  and  image 
characteristics  obtained  from  high-reso¬ 
lution  spectral  information  that  is  avail¬ 
able  only  with  HSSR  MR  imaging.  In  both 
cases,  our  results  suggest  that  HSSR  MR 


imaging  has  some  advantages  over  con¬ 
ventional  MR  imaging.  Results  of  both 
quantitative  analysis  and  evaluations 
by  experienced  radiologists  suggest  that 
HSSR  MR  images  may  show  anatomy  that 
is  more  clearly  defined  than  that  on  con- 
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Comparison  of  Fat-Suppression  Ratio,  Normalized  Surface  Area,  and  Signal-to-Noise  Ratio  in  Six  Normal  Breasts 

Fat-Suppression  Ratio* 

Normalized  Surface  Area1- 

Signal-to-Noise  Ratio 

Patient  No, 

Tl-weighted 

Fat-saturated 

MR  Images 

HSSR  MR 

Images 

T1  -weighted 
Fat-saturated 

MR  Images 

HSSR  MR 
Images 

T1  -weighted 
Fat-saturated 

MR  Images 

HSSR  MR 
Images 

1 

0,9 

3.5 

0.6 

4,9 

13 

12 

2 

1.2 

3.7 

1,4 

12,1 

31 

29 

3 

1,6 

4.2 

0.5 

1.3 

14 

30 

4 

1,5 

3.4 

2.2 

5,8 

13 

12 

5 

1,1 

3.5 

3,2 

16.6 

26 

26 

6 

Mean  ±  standard  error 

2,3 

4.0 

5,5 

5.6 

37 

27 

of  the  mean 

1 ,4  ±  0,2 

3,7  ±  0.1 

2.2  ±  0.8 

7.7  ±  2,3 

22,3  ±  4.3 

22.7  ±  3,4 

*  Fat-suppression  ratio  = 

signal  in  parenchymal  pixels  containing  primarily  water  divided  by  signal  in  pixels  containing  primarily  fat. 

T  uetinition  ot  normalized  surface  area  is  given  by  5 

in  Equation  (2).  5, 

is  designed  to  measure  the  sharpness  of  an  image. 

Figure  4,  A~Hf  Sagittal  images  of  breasts  with  carcinomas  (arrows).  Top  row:  Images  in  a 
47-year-old  patient  with  infiltrating  ductal  carcinoma,  Scarff-Bloom-Richardson  grade  II,  in  the 
left  breast.  Bottom  row:  Images  In  a  47-year-old  patient  with  infiltrating  ductal  carcinoma, 
Scarff-Bloom-Richardson  grade  II,  with  a  moderate  component  of  ductal  carcinoma  in  situ,  low 
grade,  in  the  right  breast.  A  and  E  are  T2-weighted  fast  spin-echo  MR  images  before  contrast 
medium  injection.  B  and  F  are  Tl-weighted  fat-saturated  MR  images  acquired  3  minutes  after 
contrast  medium  injection.  C  and  G  are  HSSR  MR  images  proportional  to  water  resonance  peak 
height  before  contrast  medium  injection,  D  and  H  are  water  peak-height  images  obtained  1-3 
minutes  after  the  contrast  medium  injection.  Evaluation  by  experienced  radiologists  suggested 
that  vasculature  and  lesion  edges  are  more  clearly  delineated  on  the  HSSR  MR  images  (C,  G). 


ventional  MR  images.  The  HSSR  MR  im¬ 
ages  (eg,  Fig  4)  show  well-defined  lesion 
boundaries,  which  are  important  clini¬ 
cally  because  irregular  lesion  contours  are 
a  sign  of  malignancy  (27,28).  In  addition, 
detailed  spectral  data  allowed  accurate 
fitting  of  the  water  peak  and  separation 
of  the  water  signal  from  the  spectral  base¬ 
line,  As  a  result,  small  water  signals  could 
be  accurately  detected  in  the  presence  of 
very  large  fat  signals,  even  when  the 
wings  of  fat  spectra  extended  into  the 
water  region,  owing  to  very  rapid  decay 
in  the  time  domain.  The  improved  sepa¬ 
ration  of  the  water  and  fat  signals  in¬ 
creases  sensitivity  to  small  amounts  of 


water,  and  this  increase  may  facilitate  de¬ 
tection  of  cancer. 

The  HSSR  MR  images  were  acquired 
1-3  minutes  after  contrast  media  injec¬ 
tion,  whereas  the  conventional  Tl- 
weighted  fat-saturated  MR  images  were 
acquired  3-5  minutes  after  injection, 
which  caused  some  bias  in  favor  of  HSSR 
MR  images.  However,  there  was  no  data 
acquisition  at  the  peak  of  the  contrast 
media  bolus;  both  HSSR  and  conven¬ 
tional  MR  images  were  acquired  during  a 
period  after  the  peak  enhancement  is 
usually  reached,  and  signal  intensity 
changes  during  this  period  are  generally 
slow.  Therefore,  comparison  of  the  HSSR 


and  conventional  MR  images  acquired  af¬ 
ter  contrast  media  injection  is  reason¬ 
able. 

The  HSSR  data  demonstrate  that  the 
water  line  in  breast  tissue  is  inhomoge- 
neously  broadened  and  often  contains 
multiple  resolved  components.  It  is  very 
unlikely  that  this  spectral  inhomogene¬ 
ity  is  a  consequence  of  poor  shimming. 
In  these  very  small  image  pixels,  mag¬ 
netic  field  gradients  due  to  shimming  are 
approximately  linear  and,  therefore,  broaden 
the  water  line  but  do  not  cause  detectable 
asymmetry  or  multiple  resolved  compo¬ 
nents.  Spectrally  inhomogeneous  effects 
of  the  contrast  media  were  detected  in 
many  pixels;  this  finding  is  consistent 
with  those  in  previous  studies  of  animal 
models  of  cancer  (1,3,8,22)  and  with  the¬ 
oretic  predictions  (17,19,20), 

We  sometimes  found  strong  but  op¬ 
posing  changes  In  the  amplitudes  of  dif¬ 
ferent  spectral  components  in  the  same 
image  pixel.  Thus,  high  spectral  resolu¬ 
tion  makes  it  possible  to  detect  strong 
effects  of  contrast  media  that  are  missed 
with  conventional  MR  imaging.  Large 
changes  in  a  single  small  component  of 
an  inhomogeneously  broadened  reso¬ 
nance  may  reflect  subpixel  regions  of 
high  vascular  density  associated  with  an¬ 
giogenesis.  This  hypothesis  is  consistent 
with  the  fact  that  spectrally  inhomoge¬ 
neous  contrast  media  effects  appear  to 
cluster  near  the  edges  of  tumors,  where 
strong  angiogenic  activity  is  expected 
(Fig  5,  D;  references  1  and  3). 

With  the  current  relatively  slow  method 
of  data  acquisition,  high-resolution  EPS 
MR  images  may  be  useful  for  improved 
characterization  of  suspicious  lesions 
identified  with  conventional  MR  imag¬ 
ing  or  other  modalities.  However,  much 
faster  versions  of  the  method  can  be  im¬ 
plemented  that  involve  the  use  of  two-  or 
three-dimensional  trajectories  through  k 
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Figure  5.  Effects  of  contrast  medium  on  HSSR  MR  images.  These  sagittal  images  were  obtained 
of  the  lesion  shown  in  Figure  4,  A-D.  A,  Changes  in  the  image  synthesized  from  the  first  echo  in 
the  HSSR  echo  train  are  depicted.  This  image  simulates  a  conventional  T1 -weighted  subtraction 
MR  image.  B,  Changes  in  peak  height  of  the  water  magnitude  spectra  are  depicted.  C,  Changes  in 
water  resonance  frequency  (in  hertz)  are  depicted.  Increase  in  frequency  gives  rise  to  higher  signal 
intensity  on  the  image.  D,  Pixels  are  highlighted  in  which  the  difference  spectra  (postcontrast 
spectrum  minus  precontrast  spectrum)  show  spectrally  heterogeneous  changes. 


space  during  the  proton  free  induction 
decay  (15).  Other  methods  for  more  effi¬ 
cient  and  rapid  data  acquisition  are  also 
available.  Methods  (eg,  simultaneous  ac¬ 
quisition  of  spatial  harmonics  [29])  for 
parallel  data  acquisition  can  dramatically 
increase  speed,  and  spectral  resolution 
can  be  reduced  at  the  edges  of  k  space. 
Thus,  it  is  likely  that  high-quality  HSSR 
MR  images  can  be  obtained  more  rapidly 
or  from  a  large  volume  of  tissue  with 
adequate  signal-to-noise  ratio. 

A  definitive  evaluation  of  HSSR  MR  im¬ 
aging  would  include  direct  comparison 
with  optimal  multiple-point  Dixon  meth¬ 
ods  (9-11),  optimal  fat-suppression  meth¬ 
ods  (eg,  rotating  delivery  of  excitation  off 
resonance  [30]),  and  other  MR  imaging 
methods  (eg,  magnetization  transfer  con¬ 
trast  material- enhanced  MR  imaging  [31]). 
Such  an  evaluation  is  feasible  since  EPS  MR 
imaging  or  related  methods  are  available  in 
a  number  of  laboratories  where  they  are 
used  to  improve  the  accuracy  and  sensitiv¬ 
ity  of  anatomic  and  functional  MR  imag¬ 
ing  (2,4,7,18,21,32-34).  In  addition,  the 
hardware  and  software  required  to  produce 
HSSR  EPS  data  sets  are  either  commonly 
available  or  can  be  implemented  easily. 
Thus,  more  comprehensive  clinical  testing 
of  this  approach  for  MR  imaging  of  the 
breast  is  practical. 

In  conclusion,  these  preliminary  results, 
as  well  as  findings  in  previous  work  (1,3,6- 
8,22,24,26,33),  suggest  potential  advan¬ 
tages  of  HSSR  MR  imaging  of  the  breast. 
Advantages  of  HSSR  MR  imaging  include 
improved  anatomic  accuracy  and  edge  de¬ 
lineation,  improved  fat  suppression,  in¬ 
creased  sensitivity  to  small  amounts  of  wa¬ 
ter-containing  tissue,  strong  T2*  contrast, 
increased  sensitivity  to  contrast  media,  and 
increased  sensitivity  to  subpixel  environ¬ 
ments  (1,3).  Future  improvements  in  HSSR 
MR  imaging  will  reduce  data  acquisition 
time.  Even  without  these  improvements, 
we  believe  the  information  provided  with 
HSSR  MR  imaging  may  justify  the  longer 
acquisition  times. 
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Purpose:  To  acquire  high  spectral  and  spatial  resolution 
{HISS}  MR  images  of  the  water  resonance  in  rat  brain, 
evaluate  the  lineshape  of  the  water  resonance  in  small 
voxels,  and  compare  images  derived  from  HiSS  data  with 
conventional  images. 

Materials  and  Methods:  Spectroscopic  images  of  rat  brain 
were  obtained  at  4.7  Tesla  using  phase  encoding  gradients 
only.  Spectral  resolution  in  each  voxel  was  ~8  Hz  and 
bandwidth  was  1000  Hz,  Spatial  resolution  was  —250  mi¬ 
crons  in  1-mm  slices.  Images  were  synthesized  to  show  the 
water  signal  integral,  peak  height,  linewidth,  resonance 
frequency,  and  asymmetry. 

Results:  Two  or  more  resolved  components  of  the  water 
resonance  were  detected  in  - 14%  ±  6%  of  voxels  in  the 
brains  of  eight  rats.  The  water  resonances  in  -20%  ±  10% 
of  voxels  (n  -  8}  were  highly  asymmetric.  Images  with  in¬ 
tensity  proportional  to  water  signal  peak  height,  T2*,  or  to 
selected  components  of  the  water  resonance  showed  fea¬ 
tures  that  were  not  evident  in  conventional  images. 

Conclusions:  The  complexity  of  the  water  signal  reflects 
the  anatomy  and  physiology  of  the  sub-voxelar  environ¬ 
ment,  and  maybe  a  useful  source  of  image  contrast.  HiSS 
imaging  of  brain  provides  accurate  anatomic  information, 
and  may  improve  image  contrast  and  delineation  of  subtle 
anatomic  features. 

Key  Words:  neuro-imaging:  brain  imaging;  spectroscopic 
imaging;  image  contrast;  water  resonance 
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BRAIN  IS  AN  OPTIMAL  TARGET  for  MRI  in  part  because 
the  water  resonance  in  brain  is  relatively  narrow  and 
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Lorentzian,  This  is  also  an  advantage  for  functional 
imaging  studies  where  very  small  changes  in  water  sig¬ 
nal  linewidth  must  be  detected  (1,2).  However,  there  is 
both  theoretical  (3,4)  and  experimental  (5-7)  evidence 
that  the  water  resonance  in  small  image  voxels  in  brain 
is  inhomogeneously  broadened.  An  important  source  of 
spectral  heterogeneity  is  deoxyhemoglobin,  which  pro¬ 
duces  microscopic  magnetic  susceptibility  gradients 
around  veins,  venules,  and  capillaries.  This  is  demon¬ 
strated  by  a  number  of  elegant  theoretical  studies  that 
have  modeled  BOLD  (blood  oxygen  level  dependent) 
contrast  and  changes  in  BOLD  contrast  during  stimu¬ 
lus  (e.g,,  3,4).  Experimental  work  in  this  laboratoiy  has 
shown  (8)  that  different  Fourier  components  of  water 
resonances  in  small  voxels  respond  differently  to 
changes  in  tumor  oxygenation,  Reichenbach  et  al  took 
advantage  of  magnetic  susceptibility  differences  be¬ 
tween  large  veins  and  surrounding  tissue  to  produce 
“MR  venograms”  using  one  and  two  point  spectroscopic 
imaging  (9).  Other  sources  of  anatomic  heterogeneity, 
i.e,,  gray  matter  vs,  white  matter  vs,  CSF  (cerebrospinal 
fluid)  may  produce  more  subtle  variations  in  magnetic 
susceptibility  (10-12), 

Despite  the  importance  of  MRI  for  clinical  neurology 
and  neuroscience,  and  the  dependence  of  MR  images 
on  the  structure  of  the  water  spectrum,  there  is  little 
information  concerning  the  water  signal  in  brain.  This 
paper  reports  a  preliminaiy  study  using  high  spectral 
and  spatial  resolution  to  characterize  the  water  signal 
in  very  small  voxels.  Pure  absorbance  spectra  were  cal¬ 
culated  to  improve  the  effective  spectral  resolution.  The 
results  demonstrate  that  water  resonances  in  many 
small  voxels  in  rat  brain  are  complex  with  two  or  more 
resolvable  components.  Images  synthesized  from  high 
spectral  and  spatial  resolution  (HiSS)  datasets  suggest 
that  information  inherent  in  spectral  lineshapes  pro¬ 
vides  novel  contrast  and  anatomic  detail  in  MR  images 
(11,13-16), 

MATERIALS  AND  METHODS 
Experimental  Protocol 

The  experimental  protocol  was  approved  by  the  Univer¬ 
sity  of  Chicago  Animal  Care  and  Use  Committee  (Pro¬ 
tocol  #53001),  The  brains  of  adult  Copenhagen  and 
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Fisher  rats  weighing  175-210  g  were  imaged  at  4,7 
Tesla  using  a  GE/Bruker  Omega  (Fremont,  CA)  scan¬ 
ner  equipped  with  self-shielded  gradient  coils.  During 
the  scans,  rats  were  anesthetized  by  continuous  LP. 
ketamine  (5  mg/ 100  g  per  hour)  and  rompun  (0. 1  mg/ 
100  g  per  hour).  Body  temperature  was  maintained  at 
37°C  using  a  warm  water  blanket  and  warm  air. 

An  8-element  birdcage  coil  (4.0  cm  in  diameter  and 
4.0  cm  long)  was  used  for  signal  excitation  and  detec¬ 
tion,  The  rat  head  was  secured  with  tape  in  the  center  of 
the  coil.  Localized  shimming  was  performed  on  a  5-mm 
slice  including  the  1-ram  slice  that  was  later  imaged. 
Optimal  shimming  is  important  (17)  in  HISS  imaging  so 
that  spectral  details  that  reflect  local  anatomy  and 
physiology  are  not  obscured  by  broadening.  To  mini¬ 
mize  the  effects  of  eddy  currents  on  water  signal  line- 
shape,  and  improve  separation  of  spectral  and  spatial 
information  (13,15),  high  resolution  spectroscopic  im¬ 
ages  were  acquired  using  phase  encoding  only,  with 
detection  of  the  proton  FID  (free  induction  decay)  in  the 
absence  of  applied  gradients  (18,19).  SI  (spectroscopic 
image)  data  were  acquired  using  a  1282  data  matrix, 
l~mm  slice  thickness,  and  32-mm2  field-of-view.  The 
FID  was  sampled  for  128  msec  to  obtain  spectral  reso¬ 
lution  of  7.8  Hz  and  spectral  bandwidth  of  1000  Hz  (128 
spectral  points).  Small  flip  angles  (<25°)  were  used  to 
reduce  stimulated  echoes.  The  total  acquisition  time 
required  for  each  spectroscopic  image  was  ~40  min¬ 
utes,  Data  were  processed  by  three-dimensional  Fou¬ 
rier  transform  without  apodization  to  obtain  a  high  res¬ 
olution  proton  spectrum  in  each  voxel.  Only  the  water 
resonance  was  used  to  construct  images;  minor  signals 
due  to  lipids  were  excluded  from  analysis. 
Conventional  spin  echo  (TR;  2000-4000  msec,  TE: 
10-75  msec)  and  gradient  echo  (TR:  400-2000  msec, 
TE;  10-30  msec)  images  were  acquired  with  the  same 
spatial  resolution  as  spectroscopic  images.  The  total 
acquisition  time  required  for  each  conventional  image 
was  —4  to  -12  minutes. 


Calculation  of  a  Frequency  Correction  to  Minimize 
Truncation  Errors 

Truncation  effects  often  occur  in  HiSS  datasets  because 
the  T2*  of  individual  components  of  the  proton  FID  in 
very  small  voxels  is  sometimes  long.  The  truncation 
effects  were  minimized  by  digitally  shifting  the  fre¬ 
quency  of  the  water  (or  fat)  proton  resonance  in  each 
pixel  until  the  peak  frequency  of  the  major  component 
of  the  water  signal  coincided  with  one  of  the  Fourier 
components  of  the  discrete  Fourier  transform.  At  this 
point,  the  amplitude  of  the  water  peak  in  the  Fourier 
domain  is  maximized.  This  is  taken  to  be  the  most 
accurate  estimate  of  peak  amplitude.  In  principle,  this 
procedure  works  perfectly  only  if  the  water  signal  is  a 
single  narrow  Lorentzian,  However,  we  find  that  in  the 
case  of  the  complicated  water  signals  commonly  seen  in 
tumors,  the  digital  frequency  shifting  procedure  im¬ 
proves  the  estimate  of  the  height  of  the  primaiy  water 
peak  in  each  pixel  and  significantly  reduces  truncation 
artifacts  (20), 
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Calculation  of  Pure  Absorption  and  Dispersion 
Spectra  for  Each  Pixel 

A  zeroth-order  phasing  algorithm  was  used  to  increase 
sensitivity  to  spectral  lineshapes.  If  the  phase  of  the 
acquired  spectrum  relative  to  a  pure  absorption/dis¬ 
persion  is  8,  then  its  real  (R)  and  imaginary  (i)  compo¬ 
nents  are  related  to  the  pure  absorption  (U)  and  disper¬ 
sion  (V)  modes  by  the  following  equations: 

U  —  R*  eos(8)  4  I*  sin(8)  ,  V=  —  R  •  sin(0)  4  I*  cos(0) . 

(1) 

The  integral  of  the  pure  dispersion  signal  is  zero,  i.e,, 

=  -  2 R *  sin(0)  4  2^*  cos(e)  -  0 ,  (2) 

where  the  summation  is  over  all  of  the  points  in  the 
spectra.  Eq,  [2]  can  be  solved  for  0,  i.e,, 

9  =  tan -'(Xi/2r).  (3) 

Often,  this  value  of  0  is  not  optimal  due  to  noise  and 
complexify  of  the  water  resonance.  Therefore,  the  esti¬ 
mate  for  0,  was  further  refined  by  minimizing  the  ex¬ 
pression; 


Hmin 

4_ 

i^min  4  ^max 

Hniax 

L^max 

This  holds  because  the  minimum  of  the  absorption 
mode  signal  is  zero,  and  the  sum  of  the  minimum  and 
maximum  of  the  dispersion  mode  signal  is  zero. 

Automatic  Detection  of  Resolved  Components  of 
the  Water  Resonance  and  of  Asymmetric 
Water  Peaks 

The  number  of  resolved  components  of  the  water  reso¬ 
nance  and  the  asymmetry  of  the  spectra  were  measured 
from  magnitude  spectra  because  the  frequency  and 
phase  correction  method  (described  above)  works  well 
only  for  the  primary  component  of  the  water  resonance. 
First,  all  possible  peaks  were  located,  i.e,,  those  points 
in  the  spectrum  where  the  signal  was  larger  than  at  the 
adjacent  frequencies.  Then,  signal  amplitude  was  mea¬ 
sured  as  the  difference  between  the  peak  intensity  and 
the  local  minima  between  these  candidate  peaks.  Can¬ 
didate  peaks  were  selected  as  true  peaks  if  the  differ¬ 
ence  between  the  local  maxima  and  minima  was  five 
times  larger  than  the  root  mean  square  of  the  noise  and 
if,  in  addition,  the  peaks  were  at  least  20%  of  the  height 
of  the  main  peak.  Images  were  constructed  to  show  the 
location  of  water  spectra  with  multiple  peaks. 

To  determine  the  asymmetry  of  magnitude  spectra, 
the  two  sides  of  the  spectra  on  either  side  of  the  global 
maximum  were  compared.  The  “Similarity”  of  the  two 
sides  was  defined  as: 

Similarity  (S,,S„)  =  (Sr  S,J/  Jgsf  ,  (5) 
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where  St  and  Sh  are  vectors  composed  of  the  low  field 
and  high  field  sides  of  the  magnitude  water  signal,  re¬ 
spectively,  The  two  sides  were  considered  significantly 
different  when  “Similarity”  was  less  than  0,95. 

Synthesis  of  Images  From  HiSS  Datasets 

A  large  number  of  images  can  be  synthesized  from  HiSS 
datasets,  based  on  details  of  the  water  spectrum.  The 
following  methods  were  used  to  synthesize  images  dis¬ 
cussed  here: 

I,  Peak  height  image:  the  maximum  value  of  the 
absorption  spectrum  (Umax). 

II,  Integral  image:  the  summation  of  the  absorption 
spectrum  (XL/). 

III,  T2*  image:  An  empirical  estimate  of  the  T2*  line 
width,  which  does  not  assume  a  functional  form 
of  the  water  resonance,  was  obtained  by  dividing 
the  integral  by  the  peak  height  of  each  resonance 
m/Umax  »  1/T2*), 

IV,  Maximum  derivative  image:  the  maximum  (or 
negative  minimum)  of  the  derivative  of  the  mag¬ 
nitude  spectrum  {SfmBX  or  -Sfmi n). 

V,  Second  moment  image:  the  second  moment  of 
power  spectra  (square  of  magnitude)  about  the 
point  (frequency)  where  the  first  moment  is  zero. 
The  second  moment  of  a  power  spectrum  is  the 
mean  square  of  its  distribution, 

VI,  Images  of  the  intensity  of  a  specific  component  of 
the  water  spectrum:  This  suppresses  signal  due 
to  the  primary  component  of  the  water  signal,  and 
emphasizes  sub-voxelar  anatomic  features  that 
appear  to  have  characteristic  resonance  frequen¬ 
cies  determined  by  local  anatomic  structure  and 
physiology. 

To  minimize  noise,  for  methods  (II)  and  (HI)  the  sum¬ 
mations  were  performed  only  in  a  limited  region  around 
the  peak  resonance  frequency.  Images  were  processed 
using  IDL  software  (Research  Systems,  Inc,,  Boulder, 
CO), 

RESULTS 

To  test  effects  of  shimming  and  eddy  currents  in  the 
absence  of  motion  and  physiologic  variations,  a  de¬ 
gassed  solution  of  copper  sulphate  in  a  plastic  con¬ 
tainer  was  used  as  a  phantom.  Figure  1  shows  the 
phantom  images  obtained  with:  (A)  a  gradient  echo  (TR/ 
TE  =  1000/30  msec),  (B)  the  B0  map  derived  from  HiSS 
data,  (C)  peak  height  image  obtained  from  HISS  data¬ 
sets  without  frequency  and  phase  correction,  and  (D) 
peak  height  image  obtained  from  absorption  spectra 
with  frequency  shifting.  The  spectral  and  spatial  reso¬ 
lutions  and  field  view  were  exactly  the  same  as  for  the 
rat  brain  images.  The  truncation  artifacts— seen  as  “ze¬ 
bra  stripes”  in  Fig,  1C— were  minimized  by  the  fre¬ 
quency  shifting  technique  (Fig.  ID),  The  water  signals 
are  single,  narrow  Lorentzians  in  phantoms  of  uniform 
composition,  and  no  complex  water  peaks  were  de¬ 
tected,  even  when  shimming  was  deliberately  spoiled 
using  high  order  shims  so  that  resonance  frequency 
varied  by  — 140  Hz  across  the  sample. 
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Figure  1,  Images  of  a  copper  sulphate  phantom:  (A)  Gradient 
echo  image  (TR/TE  =  1000/30  msec),  (B)  the  Bo  map  derived 
from  HiSS  data,  (C)  peak  height  image  obtained  from  HiSS 
datasets  without  frequency  and  phase  correction,  and  (D)  peak 
height  image  obtained  from  absorption  spectra  with  frequency 
shifting. 


The  water  proton  signals  in  a  significant  percentage 
of  voxels  in  rat  brain  were  inhomogeneously  broadened 
and  often  contained  multiple  resolvable  components. 
The  algorithm  described  above  detected  two  or  more 
resolved  water  peaks  in  an  average  —14%  ±  6%  of 
voxels  in  the  brains  of  eight  rats.  Water  spectra  in  an 
average  of  —20%  ±  10%  of  voxels  (n  =  8)  were  asym¬ 
metric,  Figure  2  shows  the  distribution  of  voxels  with 
multiple  peaks  (column  A)  or  asymmetric  peaks  (col¬ 
umn  C)  in  approximately  the  same  slice  in  three  differ¬ 
ent  rats  (column  B).  The  distributions  in  the  three  rats 
are  qualitatively  similar,  and  complex  water  spectra 
appear  to  be  associated  with  specific  anatomic  features. 
They  are  most  frequently  located  at  well-defined  inter¬ 
faces  between  different  types  of  tissue,  e.g.,  in  tissues 
adjacent  to  ventricles  where  there  are  dense  blood  ves¬ 
sels.  Examples  of  water  magnitude  spectra  from  two 
different  HiSS  datasets  are  also  shown  in  Figure  2.  The 
spectra  in  each  column  are  selected  from  approximately 
the  same  location,  as  Indicated  by  the  arrows.  Spectra 
from  the  same  location  have  similar  lineshapes.  To  test 
reproducibility,  the  spectrum  from  each  voxel  in  one 
dataset  was  correlated  with  the  spectrum  from  the 
same  voxel  in  a  second  dataset  acquired  immediately 
after  the  first.  The  average  correlation  coefficient  ob¬ 
tained  in  a  representative  experiment  was  0.95  over 
entire  brain. 

The  complexity  of  the  water  signal  in  many  voxels 
reflects  anatomic  and  physiologic  heterogeneity  in  the 
local  environment,  and  thus  may  be  a  useful  source  of 
image  contrast.  Therefore,  we  evaluated  a  number  of 
methods  for  synthesizing  images  based  on  analysis  of 
the  water  signal  lineshape.  Figure  3  shows  (A)  a  spin 
echo  image  (TR/TE  =  3000/50  msec),  (B)  a  T2*- 
weighted  gradient  echo  image  (TR/TE  =  1000/30 
msec),  and  (C)  a  Tr -weighted  gradient  echo  image  (TR/ 
TE  =  500/10  msec)  compared  with  (D)  water  signal 
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Figure  2,  Axial  images  from  approximately  corre¬ 
sponding  slices  in  three  different  rats.  Column  B 
shows  water  resonance  peak  height  images  de¬ 
rived  from  HiSS  data.  Column  A  highlights  the 
distribution  of  pixels  with  multiple  water  reso¬ 
nances  and  column  C  highlights  the  distribution 
of  pixels  with  asymmetric  water  resonances.  Ex¬ 
amples  of  water  signal  magnitude  spectra  from 
two  different  HiSS  datasets  are  also  shown;  the 
locations  of  the  spectra  are  indicted  by  the  arrows. 


peak  height,  (E)  T2*,  and  (F)  integral  images — all  from 
the  same  rat  brain.  The  acquisition  times  required  for 
spin  echo  and  gradient  echo  images  were  —12,8  and 
-6,4  min,  respectively,  while  the  acquisition  time  for 
the  spectroscopic  image  was  -40  minutes.  The  peak 
height  and  T2*  images  synthesized  from  pure  absorp¬ 
tion  spectra  in  the  HiSS  datasets  show  some  features 
that  are  not  evident  in  the  conventional  images.  The 
C3F,  cerebellum,  hippocampus,  caudate  putamen,  and 
the  central  nucleus  of  the  inferior  colliculus  are  partic¬ 
ularly  clearly  defined.  The  peak  height  and  T2*  HiSS 
images  in  Figures  3D  and  3E  show  dark  lines  that  may 
correspond  to  veins,  for  example,  indicated  by  the  ar¬ 
rows,  These  features  are  not  seen  in  conventional  im¬ 
ages. 

Further  comparisons  between  conventional  images 
and  HiSS  images  are  shown  in  Figure  4  for  three  rat 
brains.  Gradient  echo  T2*~weighted  images  (first  col¬ 
umn  A)  are  compared  with  the  corresponding  peak 


height  images  from  HiSS  datasets  (second  column  B). 
In  addition,  other  images  obtained  from  HiSS  data  us¬ 
ing  methods  IV  (“maximum  derivative”),  V  ("second  mo¬ 
ment”),  and  VI  (“specific  frequency  component”),  as  de¬ 
fined  in  the  Methods  section,  are  shown  in  column  C. 
The  gradient  echo  imaging  parameters  for  rows  1  to  3 
were:  TR/TRE  =  1000/25  msec,  1000/20  msec,  and 
2000/30  msec,  respectively.  The  corresponding  acqui¬ 
sition  times  were  6,4  minutes,  4,3  minutes,  and  12,8 
minutes.  The  peak  height  images  (column  B)  provide 
improved  contrast,  SNR,  and  enhanced  detail,  e.g„  in 
the  basal  ganglion — particularly  in  row  2,  column  B, 
compared  to  T2*~weighted  gradient  echo  images  (col¬ 
umn  A),  The  maximum  derivative  image  (row  1,  column 
C)  shows  strong  contrast  where  water  spectra  have  a 
steeper  slope,  for  example,  near  the  3rd  ventricle  indi¬ 
cated  by  the  arrow.  The  second  moment  image  (row  2, 
column  C)  shows  features— particularly  near  the  ante¬ 
rior  pretectal  nucleus  indicated  by  the  arrow — that  can- 
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Figure  3.  Conventional  (A)  spin  echo,  (B)  gradi¬ 
ent  echo  T2*-weighted,  and  (C)  gradient  echo  Tr 
weighted  images  compared  with  (□)  peak  height, 
(E)  T2*,  (F)  integral  images  extracted  from  the  ab¬ 
sorption  spectra  calculated  from  HISS  datasets. 


not  be  seen  in  peak  height  and  gradient  echo  images. 
Finally,  the  image  from  a  selected  off-peak  Fourier  com¬ 
ponent  (row  3,  column  C)  shows  structures  near  den¬ 
tate  gyrus  and  hippocampal  area  indicated  by  the  arrow 
that  are  not  seen  in  gradient  echo  and  peak  height 
images. 

DISCUSSION 

These  results  show  that  many  water  resonances  in 
small  voxels  in  rat  brains  are  highly  non-Lorentzian, 
with  multiple  resolved  components.  The  non-Lorent- 
zian  spectral  features  were  identified  using  conserva¬ 
tive  criteria  that  probably  underestimate  the  prevalence 
of  complex  water  spectra  in  rat  brain.  The  results  ob¬ 


tained  in  phantoms  suggest  that  most  of  these  features 
are  not  due  to  poor  shimming  or  eddy  currents.  At  the 
spatial  and  spectral  resolution  of  the  present  experi¬ 
ments,  water  resonances  from  single  voxels  inside  uni¬ 
form  phantoms  were  single,  narrow  Lorentzian  lines. 
Poor  shimming  caused  broadening  but  multiple  and 
asymmetric  peaks  were  not  observed.  Eddy  currents 
due  to  phase  encoding  and  slice  selection  gradients 
would  have  been  roughly  similar  in  the  phantoms  and 
in  rat  brain,  but  there  were  no  observable  effects  of  eddy 
currents  in  phantoms.  This  is  consistent  with  the  ex¬ 
pectation  that  macroscopic  magnetic  field  gradients, 
e.g.,  due  to  poor  shimming  or  sinuses,  are  roughly  lin¬ 
ear  across  small  voxels,  and  thus  will  only  cause  broad- 
ening.  In  fact,  the  spectral  characteristics  of  water  res- 


Figure  4.  Comparison  of  gradient  echo  T2- 
weighted  images  (column  A)  with  peak  height  im¬ 
ages  (column  B)  and  other  images  (column  C) 
extracted  from  the  HtSS  datasets  for  three  differ¬ 
ent  rats.  In  column  C,  rows  1,  2,  and  3  are  the 
maximum  derivative”,  “second  moment",  and 
“specific  frequency  component"  images,  respec¬ 
tively. 


onances  in  brain  appear  to  be  associated  with  specific 
anatomic  features,  rather  than  with  large  macroscopic 
B0  gradients.  This  suggests  that  they  are  caused  pri¬ 
marily  by  local  variations  (i.e.,  subvoxelar  or  between 
neighboring  voxels)  in  magnetic  susceptibility  that  re¬ 
flect  the  heterogeneity  of  local  anatomy  and  physiology. 
For  example,  small  blood  vessels  or  sharp  transitions  in 
tissue  composition  (12)  cause  shifts  in  resonance  fre¬ 
quency  that  may  produce  distinguishable  components 
of  the  water  resonance  ,  Spectral  resolution  higher  than 
that  of  the  present  data  (~8  Hz)  might  show  more  detail 
than  is  present  in  the  current  images,  albeit  at  the  cost 
of  decreased  SNR.  Ultimately,  macroscopic  field  gradi¬ 
ents  will  limit  our  ability  to  detect  local  gradients  in 
magnetic  susceptibility — but  this  problem  can  be  offset 
by  use  of  sophisticated  shimming  methods  (21), 

This  report  describes  preliminary  work,  A  great  deal 
remains  to  be  done  to  identify  the  sources  of  inhomog- 
enous  broadening  in  each  region,  and  to  determine 
whether  spectra!  information  can  be  analyzed  to  produce 
novel,  clinically  useful  images.  Nevertheless,  the  data 
suggest  potential  advantages  of  high  spectral  and  spatial 
resolution  imaging  that  merit  further  examination,  HISS 
images  of  normal  rat  brain  show  features  that  are  not 
detectable  in  conventional  images.  Sensitivity  to  specific 
spectral  components  may  provide  information  at  the  sub¬ 
voxelar,  or  perhaps  microscopic  level  (8,15,22). 

Qualitative  differences  between  HISS  images  and 
conventional  images  are  probably  not  due  to  differences 
in  the  slgnal-to-noise  ratio  of  the  images.  The  acquisi¬ 
tion  times  of  HiSS  images  were  substantially  longer 
than  for  conventional  images  due  to  the  requirements 
of  the  phase  encoding  method.  However,  the  HiSS  data¬ 
sets  were  predominantly  composed  of  very  low  signal 
components  of  the  proton  FID  at  long  evolution  times 
(out  to  128  msec),  while  the  T2* -weighting  of  the  con¬ 
ventional  gradient  echo  images  was  much  less  pro¬ 
nounced  (echo  times  of  10-30  msec).  Thus,  the  overall 
signal-to-noise  ratio  of  HiSS  and  conventional  images 
was  similar.  It  is  likely  that  the  differences  between 
HiSS  and  conventional  images  derive  primarily  from  the 
fact  that  the  HiSS  images  are  much  more  heavily  T2*~ 
weighted  than  conventional  gradient  echo  images, 
without  the  distortions  that  generally  accompany 
strong  T2*-  weighting,  and  that  blurring  due  to  local 
susceptibility  gradients  can  be  minimized  in  HiSS  im¬ 
ages,  In  addition,  HiSS  imaging  can  sometimes  distin¬ 
guish  signals  from  different  anatomic  structures  based 
on  resonance  frequency. 

Echo  planar  spectroscopic  imaging  techniques  allow 
acquisition  of  HiSS  data  with  dramatically  reduced  run 
times  and  can  be  used  for  high  resolution  anatomic  and 
functional  imaging  (11,13,14,23)  in  a  clinical  setting. 
Therefore,  HiSS  MM  of  brain  can  easily  be  evaluated  in 
patients  (2), 
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